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Abstract. We prove well-posedness for the 3-D compressible Euler equations with mov- 
ing physical vacuum boundary, with an equation of state given by p{p) = C-yp'^ for 
7 > 1. The physical vacuum singularity requires the sound speed c to go to zero as the 
square-root of the distance to the moving boundary, and thus creates a degenerate and 
characteristic hyperbolic free- boundary system wherein the density vanishes on the free- 
boundary, the uniform Kreiss— Lopatinskii condition is violated, and manifest derivative 
loss ensues. Nevertheless, we are able to establish the existence of unique solutions to 
this system on a short time-interval, which are smooth (in Sobolev spaces) all the way 
to the moving boundary, and our estimates have no derivative loss with respect to initial 
data. Our proof is founded on an approximation of the Euler equations by a degenerate 
parabolic regularization obtained from a specific choice of a degenerate artificial viscosity 
term, chosen to preserve as much of the geometric structure of the Euler equations as 
possible. We first construct solutions to this degenerate parabolic regularization using 
a new higher-order Hardy-type inequality; we then establish estimates for solutions to 
this degenerate parabolic system which are independent of the artificial viscosity pa- 
rameter. Solutions to the compressible Euler equations are found in the limit as the 
artificial viscosity tends to zero. Our regular solutions can be viewed as degenerate vis- 
cosity solutions. Out methodology can be applied to many other systems of degenerate 
and characteristic hyperbolic systems of conservation laws. 
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1. Introduction 

1.1. The compressible Euler equations in Eulerian variables. For < t <T, the evo- 
lution of a three-dimensional compressible gas moving inside of a dynamic vacuum boundary 
is modeled by the one-phase compressible Euler equations: 

p[ut + u ■ Du] + Dp{p) = in n{t) , (1.1a) 

pt + div{pu) = in n{t) , (1.1b) 

p = on r{t), (1.1c) 

Vir{t)) = u-n{t) (l.ld) 

{p,u) = {po,uo) on f7(0), (l.le) 

n{o) = n. (i.if) 

The open, bounded subset ^l{t) C M'^ denotes the changing volume occupied by the gas, 
T{t) := dfl{t) denotes the moving vacuum boundary, V{T{t)) denotes the normal velocity 
of T{t), and n{t) denotes the exterior unit normal vector to T{t). The vector-field u — 
(ui, M2, U3) denotes the Eulerian velocity field, p denotes the pressure function, and p denotes 
the density of the gas. The equation of state p{p) is given by 

p{x,t) = C-y p{x,ty for 7>1, (1.2) 

where C-y is the adiabatic constant which we set to unity, and 

p > in n{t) and p on r{t) . 

Equation p. lap is the conservation of momentum; (jl.lbp is the conservation of mass; the 
boundary condition (jl.lcp states that the pressure (and hence the density function) vanish 
along the moving vacuum boundary r(i); (jl.ldp states that the vacuum boundary T{t) is 
moving with speed equal to the normal component of the fluid velocity, and (|l.lep - (jl.lf| 
are the initial conditions for the density, velocity, and domain. Using the equation of state 
pr^ . (jl.lap is written as 

p[ut + u- Du]+ Dp'' = in n{t). ( [TTal ) 



1.2. Physical vacuum. With the sound speed given by c :— ^ dp/ dp and N denoting the 
outward unit normal to the initial surface F, satisfaction of the condition 

§ < on F (1.3) 

defines a physical vacuum boundary (see [31], [53], [5S], [53], [27], [H]), where cq = c\t=o 
denotes the initial sound speed of the gas. 

The physical vacuum condition (|1.3p is equivalent to the requirement that 

^<OonF, (1.4) 

a condition necessary for the gas particles on the boundary to accelerate. Since p^ > Q in 
f2, (jl.4p implies that for some positive constant C and x & U, near the vacuum boundary F, 

pl^^{x)>Cdisi{x,T). (1.5) 

Because of condition (jl.Sp , the compressible Euler system (11.11) is a degenerate and char- 
acteristic hyperbolic system which violates the uniform Kreiss-Lopatinskii condition |18] . 
and to which standard methods of symmetric hyperbolic conservation laws are apparently 
ineffective. The moving boundary is characteristic because of the evolution law (jl.ldp . and 
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the system of conservation laws is degenerate because of the appearance of the density 
function as a coefficient in the nonhnear wave equation which governs the dynamics of the 
divergence of the velocity of the gas; in turn, weighted estimates show that this wave equa- 
tion indeed loses derivatives with respect to the uniformly hyperbolic non-degenerate case 
of a compressible liquid, wherein the density takes the value of a strictly positive constant 
on the moving boundary [7] . We provide a brief history of results in this area in Section 11.91 
below. 

We note that with a faster rate of degeneracy of the density function, such as, for example, 
dist(x, T{t))'' for 6 = 2,3, the analysis becomes significantly easier; for instance, if 6 = 2, 

Dp'^~^ (x t) 

then J _i is bounded for all x G fl. This bound makes it possible to readily control 

error terms in energy estimates, and in effect removes the singular behavior associated 
with the physical vacuum condition (jl.5l) . On the other hand, if p^~^ tends to zero like 
dist(x, r(i))'' for 6 = 2,3, then the gas cannot accelerate into vacuum. 

1.3. Fixing the domain and the Lagrangian variables on fi. We transform the system 
()1.1|) into Lagrangian variables. We let ri{x, t) denote the "position" of the gas particle x at 
time t. Thus, 

dti] — uo ri for t > and 'q{x, Q) = x 
where o denotes composition so that [u o ri]{x, t) :— u{ri{x, t),t) . We set 
V — uo rj (Lagrangian velocity), 
f = por] (Lagrangian density) , 
A = [D'r]]~^ (inverse of deformation tensor), 
J = det Dt] (Jacobian determinant), 
a — J A (transpose of cofactor matrix). 



Using Einstein's summation convention defined in Section [2. SI below. and using the notation 
F,k to denote -ir-, the ki/i-partial derivative of F for k — 1,2, 3, the Lagrangian version of 
equations (|l.la|) - (jl.lbp can be written on the fixed reference domain f2 as 

K + Afr,fe=0 inr!x(0,T], (1.6a) 

ft + fAlv\,^0 mnx{0,T], (1.6b) 

/ = mnx{0,T], (1.6c) 

{f,v,r]) ^ {po,uo,e) in x {t = 0} , (1.6d) 

where e{x) — x denotes the identity map on fl. 

Since Jt = JAjv^j and since J(0) = 1 (since we have taken r]{x, 0) = x), it follows that 

f = PoJ-\ (1.7) 

so that the initial density function po can be viewed as a parameter in the Euler equations. 
Let r := dfl denote the initial vacuum boundary; using, that A^ — a^, we write the 
compressible Euler equations (II. 6p as 

povl + a1{plJ~^),k = in r! X (0,T] , (1.8a) 

{ri,v) ^ {e,Uf)) \nVLx{t = 0}, (1.8b) 

pT'^ = Q onT, (1.8c) 
with pg ^(x) > C dist(a;, F) for x G near F. 
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1.4. Setting 7 = 2. We will begin our analysis for the case that 7 = 2, and in Section [Tl] 
we will explain the modifications required for the case of general 7 > 1. 

With 7 set to 2, we thus seek solutions 77 (t) to the following system: 

povl + a'l{plJ-^),k =0 in X (0,T] , (1.9a) 

(r/,w) = (e,ito) onr2x{t = 0}, (1.9b) 
po = on r , (1.9c) 

with po{x) > C dist(a;, F) for x d Q near T. 
The equation (|1.9a|) is equivalent to 

vl + 2A'lipoJ-'),k^0, (1.10) 

and p.lOp can be written as 

vl + Poa'^J-^k +2po,k a'lJ-^ = . (1.11) 
Because of the degeneracy caused by po = on F, all three equivalent forms of the com- 
pressible Euler equations are crucially used in our analysis. The equation (jl.9ap is used for 
energy estimates, while (ll.lOp is used for estimates of the vorticity, and (|l.lip is used for 
additional elliptic-type estimates used to recover the bounds for normal derivatives. 

1.5. The reference domain Q.. To avoid the use of local coordinate charts necessary for 
arbitrary geometries, and to simplify our exposition, we will assume that the initial domain 

C M'^ at time t = is given by 

n = {{xuX2,xz)^W' I [xi,X2)(.T'', x^e (0,1)}, 

where denotes the 2-torus, which can be thought of as the unit square with periodic 
boundary conditions. This permits the use of one global Cartesian coordinate system. At 
t = 0, the reference vacuum boundary is comprised of the bottom and top of the domain f2 
so that 

r = {X3 = 0} U {X3 - 1} . 
Then, according to the evolution law for the moving vacuum boundary T(t) given by 
(jl.ldp . we have that 

(We will sometimes write ri{t, F) to denote ri{t){r).) Hence, solving (|1.9I) for ri{t) (and v{t) = 
Tit{t)) completely determines the motion and regularity of the moving vacumn bomidary T(t). 

1.6. The higher-order energy function for the case 7 = 2. The physical energy 
/n[^'^ol^l^ Po'^~^]dx is a conserved quantity, but is far too weak for the purposes of 
constructing solutions; instead, we consider the higher-order energy function 

4 

+ II curl,, v{t)\\l + WpoB^ curl,, v{t)\\l , (1.12) 

where B = ^gf^i sfj) curl,, v — [curlu] o rj. Section [5] explains the notation. 

We also define Mq = P{E{0)), where P denotes a polynomial function of its argument. 

While the higher-order energy function E{t) is not conserved, we will construct solutions 
to (|1.9|) for which supjgjQ -B(t) remains bounded whenever T > is taken sufficiently 
small; the bound depends only on E{0). 
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1.7. Main Result. 

Theorem 1.1 (Existence and uniqueness for the case 7 = 2). Suppose that pq e H'^{^1), 
Po{x) > for X d fl, and po satisfies il.5\} . Furthermore, suppose that uq is given such that 
Mq < 00. Then there exists a solution to m.9\) (and hence lil.l]) ) on [0,T] for T > taken 
sufficiently small, such that 

sup Bit) < 2Mq . 
te[o,T] 

In particular, the flow map 77 G L°"{0,T; H*{fl)) and the moving vacuum boundary T{t) is 
of Sobolev class H^'^. 

Moreover if the initial data satisfies 

a=0 

+ \\cni\^v{0)\\l + \\pQd^cnTl^v{0)\\l < 00 , (1.13) 
then the solution is unique. 

Remark 1.2. The case of arbitrary 7 > 1 is treated in Theorem l 1 1 . 1\ below. 

Theorem 1 1 . 1 1 also covers the 2-D case that C K^. We established the analogous result 
in 1-D in [11]. We note that by using a collection of local coordinate charts, we could 
modify our proof to allow for arbitrary initial domains J7, as long as the initial boundary is 
of Sobolev class H^-^. 

The multi-D physical vacuum problem is not only a characteristic hyperbolic system, 
but it is also degenerate because the density function vanishes on the boundary F. In 1- 
D, the two characteristic curves of the isentropic system intersect with the moving vacuum 
boundary T{t) tangentially; this triple point of intersection is suggestive of singular behavior. 
While the degeneracy produces "honest" derivative loss with respect to uniformly hyperbolic 
systems, we develop a methodology based on nonlinear estimates which provides us with a 
priori control of smooth solutions which do not suffer from the derivative loss phenomenon 
(see [Hj for the a priori estimates to this problem). As we will outline below, our method for 
constructing smooth solutions does not rely on linearization, Kreiss-Lopatinskii theory, or 
the Nash-Moser iteration scheme, but rather on a carefully chosen nonlinear approximation 
to the characteristic and degenerate Euler equations, which preserves a great deal of the 
nonlinear structure of the original system. 

1.8. History of prior results on the analysis of multi-D free-boundary Euler prob- 
lems. 

1.8.1. The incompressible setting. There has been a recent explosion of interest in the anal- 
ysis of the free-boundary incompressible Euler equations, particularly in irrotational form, 
that has produced a number of different methodologies for obtaining a priori estimates, 
and the accompanying existence theories have mostly relied on the Nash-Moser iteration to 
deal with derivative loss in linearized equations when arbitrary domains are considered, or 
complex analysis tools for the irrotational problem with infinite depth. We refer the reader 
to [2], [10], [20], [22], [33], [35], [41], [42], [43], and [46] for a partial list of papers on this 
topic. 
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1.8.2. The compressible setting. The mathematical analysis ol moving hypersurfaces in the 
multi-D compressible Euler equations is essential for the understanding of shock waves, 
vortex sheets or contact discontinuities, as well as phase transitions such as the motion of 
gas into the vacuum state considered herein. 

The stability and regularity of the multi-D shock solution was initiated in [29] and ex- 
tensively studied by [13], [14], [15], and [32] (see the references in these articles for a more 
extensive bibliography). The shock wave problem is non-characteristic on the boundary, 
and in fact, produces the so-called dissipative boundary conditions, and satisfies the uni- 
form Kreiss-Lopatinskii condition. Even so the methodologies employed produce derivative 
loss with respect to initial data. 

More delicate than the non-characteristic case, is the characteristic boundary case, en- 
countered in the study of vortex sheet or current vortex sheet problems. This class of 
problems has been studied by [3, [S], [6], [38], [40] and others, and has the relative dis- 
advantage of violating the uniform Kreiss-Lopatinskii condition, which produces derivative 
loss in the linearization, similar to that experienced by many authors in the incompressible 
flow setting (both irrotational flows and flows with vorticity). 

1.9. History of prior results for the compressible Euler equations with vacuum 
boundary. The physical vacuum free-boundary problem, also described as the physical 
vacuum singularity, has a rich history, as well as a great deal of renewed interest (see [I]). 

Some of the early developments in the theory of vacuum states for compressible gas dy- 
namics can be found in j27j and |21j . We are aware of only a handful of previous theorems 
pertaining to the existence of solutions to the compressible and undamped Euler equation|3 
with a moving vacuum boundary. In j30j . compactly supported initial data was considered, 
and the compressible Euler equations were treated as a PDE set on x (0,T]. Unfortu- 
nately, with the methodology of [3^, it is not possible to track the location of the vacuum 
boundary (nor is it necessary); nevertheless, an existence theory was developed in this con- 
text, by a variable change that permitted the standard theory of symmetric hyperbolic 
systems to be employed, but, the constraints on the data were too severe to allow for the 
evolution of the physical vacuum boundary. 

Existence and uniqueness for the 3-D compressible Euler equations modeling a liquid 
rather than a gas was established in [23]. As discussed in [7], for a compressible liquid, the 
density p > A > is assumed to be a strictly positive constant on the moving vacuum bound- 
ary r(i) and p is thus uniformly bounded from below by a positive constant. As such, the 
compressible liquid provides a uniformly hyperbolic, but characteristic, system. Lagrangian 
variables combined with Nash-Moser iteration was used in [53] to construct solutions. More 
recently, [39] provided an alternative proof for the existence of a compressible liquid, employ- 
ing a solution strategy based on symmetric hyperbolic systems combined with Nash-Moser 
iteration, but as stated in Remark 2.2 of that paper, the 7-gas law equation-of-state p — 
cannot be used. 

In the presence of damping, and with mild singularity, some existence results of smooth 
solutions are available, based on the adaptation of the theory of symmetric hyperbolic 
systems. In [25], a local existence theory was developed for the case that c" (with < a < 1) 
is smooth across F, using methods that are not applicable to the local existence theory for 
the physical vacuum boundary. An existence theory for the small perturbation of a planar 



The parabolic free-boundary viscous Navier-Stokes equations do not experience the same sort of analyt- 
ical difficulties as the compressible Euler equations, so we do not focus on the viscous regime in this paper. 
We refer the reader to |16| . | 28| . |31| , and |34| for the analysis of the corresponding viscous system. 
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wave was developed in [44]. See also j26] and [45j, for other features of the vacuum state 
problem. 

In the 1-D setting, recently, [T7] have established existence and uniqueness using weighted 
Sobolev norms for their energy estimates. From these weighted norms, the regularity of the 
solutions cannot be directly determined. Letting d denote the distance function to the 
boundary dl, and letting || • ||o denote the L^(r2)-norm, an example of the type of bound 
that is proved for their rescaled velocity field u in 17 is the following: 



This bound is obtained from their paper by considering the case 7 = 3, fc = 1, and making 
the assumption that (j) — ^ (using the variable terminology of their paper) which is certainly 
true near the boundary. The problem with inferring the regularity of u from this bound can 
already be seen at the level of an H^{i}) estimate. In particular, the bound on the norm 
+ 2^2: II Q only implies a bound on lldu^rcHo and Hu^Ho if the integration by parts on 
the cross-term, 



can be justified, which in turn requires having better regularity for than the a priori 
bounds provide. Any methodology which seeks regularity in (unweighted) Sobolev spaces 
for solutions must contend with this type of issue. 

We overcame this difficulty in 1-D [TT] by constructing (sufficiently) smooth solutions to 
a degenerate parabolic regularization and consequently avoiding this sort of integration-by- 
parts difficulty. Our solution strategy in was based on a 1-D version of our higher-order 
Hardy-type inequality. In this paper, we extend our ideas to the multi-D setting. 

1.10. Outline of the paper and our methodology. Section [2] defines the notation used 
throughout the paper. In Section [31 we introduce our higher-order Hardy inequality for 
functions on J7 that vanish on F; this inequality is of fundamental importance to our strategy 
for constructing solutions. In this section, we also state a lemma on K-independent estimates 
for equations Kft + f = g, which will be of great use to us in the elliptic-type estimates 
that we shall employ for bounding normal derivatives. We end this section with a standard 
weighted embedding into standard Sobolev spaces. Section 2] defines the Lagrangian curl and 
divergence operators. In Section [SI we provide basic differentiation rules for the Jacobian 
determinant J and cofactor matrix a, and state the basic geometric and relevant analytical 
properties of the cofactor matrix. Section [5] provides some well-known elliptic estimates 
based on the Hodge decomposition of vector fields, as well as some basic trace estimates for 
the normal and tangential components of vectors fields in L^(ri). 

In Section [71 we introduce the degenerate parabolic approximation (|7.2p to the compress- 
ible Euler equations (|1.9p . which takes the form povl + a^{plJ~^),k +ndt [a*^(/9§ J~^),fc ] ~ 0, 
where k > denotes the artificial viscosity parameter, and with the special choice of the de- 
generate parabolic operator Kdt [af^(po J~^),fe ] which preserves a majority of the geometric 
structure of the Euler equations. In particular, the structure of the energy estimates for the 
horizontal space derivatives as well as time derivatives are essentially preserved, the elliptic- 
type estimates for vertical (or normal) derivatives are kept intact, while the estimates for 
vorticity are not exactly preserved, but can still be obtained with some additional structural 
observations employed. 

Section[51is devoted to the construction of solutions to the degenerate parabolic K-problcm 
(j7.2p on a time interval [0,Tk], where may a priori approach zero as k — )• 0. The 1-D 



\\du\\l + \\dux\\l + \\duxx + 2Ma;||o + \\du.xxx + "^-Uxx - 2d ^ u^Wl 

+ \\dUxxxx + -iUxxx - 4rf"^ UxxWl < 00 • 



(1.14) 
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version of the parabolic K-problem has been studied by us in 11^ and also in [12 in the 
context of the Wright-Fisher diffusion arising in mathematical biology. 

The construction of solutions in the 3-D setting is significantly more challenging. Our 
approach is to (1) compute the Lagrangian divergence of the K-problem to find a nonlin- 
ear degenerate parabolic equation for podiv^t;, (2) compute the Lagrangian curl of the 
K-problem to find the evolution equation for curl^ v, and (3) to consider the vertical (or 
normal) component of the trace of the K-problem on the boundary T, and find an evolution 
equation for i;"^. We then linearize these three evolution equations, and obtain a solution to 
the linearized problem via an additional approximation scheme, which requires us to hori- 
zontally smooth the linearized boundary evolution PDE for , using convolution operators 
on r. We find a fixed-point to this horizontally smoothed problem using the contraction 
mapping principle, and then perform energy estimates to find a solution on a time-interval 
which is independent of the horizontal convolution parameter. An additional contraction 
mapping argument is then made to find a solution of the nonlinear K-problem. One of 
the serious subtleties of our analysis involves the solution and regularity of the degenerate 
parabolic equation for the po div^ v. 

In Section [9l we establish K-independent estimates for the solutions that we have con- 
structed to the K-problem (|7.2|) . This is done by a combination of energy estimates for the 
horizontal and time-derivatives of ?7(i), which rely on the determinant structure of the Euler 
equations in Lagrangian variables, followed by elliptic-type estimates that give bounds on 
the vertical derivatives of r]{t) and its time-derivatives. 

Section [TOl uses these k- independent estimates to construct a solution to the compressible 
Euler equations as a limit of the sequence of parabolic solutions as k — > 0. Uniqueness is 
proven as well. 

Finally, in Section 1111 we describe the modifications which are necessary for the case of 
general 7 > 1. 

The methodology developed for the multi-D compressible Euler equations with physical 
vacuum singularity is somewhat general, and can be applied to a host of other degenerate 
and characteristic hyperbolic systems of conservation laws such as the equations of magneto- 
hydrodynamics. 



1.11. Generalization of the isentropic gas assumption. The general form of the com- 
pressible Euler equations in three space dimensions are the 5x5 system of conservation 
laws 

p[ut + u ■ Du] + Dp{p) = , (1.15a) 
pt + div{pu) ^ , (1.15b) 
{pe)t+dW{pu€ + pu)=0, (1.15c) 

where (jl.lSal) . (jl.l5bp and (jl.lScI) represent the respective conservation of momentum, mass, 
and total energy. Here, the quantity € is the sum of contributions from the kinetic energy 
^litp, and the internal energy e, i.e.,£ = jjup -I- e. For a single phase of compressible liquid 
or gas, e becomes a well-defined function of p and p through the theory of thermodynamics, 
e = e{p,p). Other interesting and useful physical quantities, the temperature T{p,p) and 
the entropy S{p,p) are defined through the following consequence of the second law of 
thermodynamics 

T dS ~ de — ^ dp . 

p^ 
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For ideal gases, the quantities e, T, S have the exphcit formulae: 

p T 



Pil - 1) 7 - 



T{p,p) = ^ 



p — e , 7 > 1, constant . 

In regions of smoothness, one often uses velocity and a convenient choice of two additional 
variables among the five quantities S, T,p, p, e as independent variables. For the Lagrangian 
formulation, the entropy S plays an important role, as it satisfies the transport equation 

St + iu- D)S = 0, 

and as such, S o ri — Sq, where Sq{x) = S{x,0) is the initial entropy function. Thus, by 
replacing / with e^°^ p^J''' , our analysis for the isentropic case naturally generalizes to the 
5x5 system of conservation laws. 

2. Notation and Weighted Spaces 

2.1. The gradient and the horizontal derivative. The reference domain J7 is defined 
in Section 11.51 Throughout the paper the symbol D will be used to denote the three- 
dimensional gradient vector 

Odd 
dxi ' dx2 ' dxj 
and let d denote the horizontal derivative 

d_ _d_ 

dxi ' dx2 



D 
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2.2. Notation for partial differentiation. The kih partial derivative of F will be denoted 

2.3. The divergence and curl operators. We use the notation divV^ for the divergence 
of a vector field V on fi: 

di-vV ^V\i+V^,2+V\3 , 
and we use curly to denote the curl of a vector V on fl: 

curly = (y'\2-y^3 ,v\3-v\i ,v\i-v\2) ■ 

Throughout the paper, we will make use of the permutation symbol 

{1, even permutation of {1, 2, 3}, 
— 1, odd permutation of {1, 2, 3}, (2.1) 
0, otherwise , 

This allows us to write the ith component of the curl of a vector-field V as 

[curiy]i — £ijkV^ ,j or equivalently curlF — e.jkV'',j 

which agrees with our definition above, but is notationally convenient. 

We will also define the Lagrangian divergence and curl operators as follows: 

div„W^AlW\j , (2.2) 

and 

cml^V = e.,kA^jW\r . (2.3) 
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In the sequel we shall also use the notation divfj and curlf^ to mean the operations defined 
by (|2.2p and (|2.3|) . respectively, with A = [Dtj]^^ replacing A, 

Finally, we will make use of the 2-D divergence operator divr for vector-fields F on the 
2-D boundary F: 

div^F = F^l+i^^2 . (2.4) 



2.4. Sobolev spaces on f2. For integers fc > and a smooth, open domain ft of R'^, we 
define the Sobolev space H^{n) R^)) to be the completion of C°°{n) (C°°(rj;R3)) 

in the norm 

\\^h--= I E / \D''u{x)\'dx 

for a multi-index a G Z'^, with the standard convention that \a\ ~ 01-1-02-1-03. For real 
numbers s > 0, the Sobolev spaces H^{il.) and the norms || • \\s are defined by interpolation. 
We will write H''{il) instead of i7^(f2;R'^) for vector- valued functions. In the case that 
s > 3, the above definition also holds for domains $7 of class H^. 

Our analysis will often make use of the following subspace of H^{V,): 

Hq{V,) = {u & H^{i}) : w = on F , (xi, X2) H' X2) is periodic } , 

where, as usual, the vanishing of w on F is understood in the sense of trace. 

We will on occasion also refer to the Banach space W^'°° (fi) consisting of L°° (17) functions 
whose weak derivatives are also in L°°{il). 




2.5. Einstein's summation convention. Repeated Latin indices i, j, fc,, etc., are summed 
from 1 to 3, and repeated Greek indices a, /3, 7, etc., are summed from 1 to 2. For example, 

r^i Talir'i . v^2 ^^2 3F' jaB 3G' 



Pni ■- I]i=l,3 dx.dxi ' ^^'i I°''^G\p '— J2i=l J2a=l Y^jS^ 



1 dXa dXj3 ' 



2.6. Sobolev spaces on F. For functions u e H''{T), fc > 0, we set 




for a multi-index a G Z^. For real s > 0, the Hilbert space iJ*(F) and the boundary norm 
I • Is is defined by interpolation. The negative-order Sobolev spaces iJ^*(r) are defined via 
duality: for real s > 0, 
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2.7. Notation for derivatives and norms. Throughout the paper, we will use the fol- 
lowing notation: 

^ , . , f d d d 

D — three-dimensional gradient vector = — — , — — , — — 

\dxi 0x2 0x3 

f d d 

d = two-dimensional gradient vector or horizontal derivative = — — , — — 

\oxi 0x2 

div = three-dimensional divergence operator , 
div^ = three-dimensional Lagrangian divergence operator , 
curl = three-dimensional curl operator , 
curl,; = three-dimensional Lagrangian curl operator , 
divr = two-dimensional divergence operator , 
II • lis = H'^{Q) interior norm, 

I • Is = H^{r) boundary norm. 

2.8. The outward unit normal to T. We set N = (0, 0, 1) on {xg = 1} and N = (0, 0, -1) 

on {x3 — 0}. We use the standard basis on W": e\ = (1, 0, 0), 62 = (0, 1, 0) and 63 = (0, 0, 1). 

3. A HIGHER-ORDER HARDY-TYPE INEQUALITY AND SOME USEFUL LEMMAS 

We will make fundamental use of the following generalization of the well-known Hardy 
inequality to higher-order derivatives: 

Lemma 3.1 (Higher-order Hardy- type inequality). Let s > I be a given integer, and suppose 
that 

If d{x) > for X € fl, d E H^{VL), r — max(s — 1, 3), and d is the distance function to dVl 

near dil, then — G iJ^^^ffi) and 
d 



<C\\u\\,. (3.1) 

s-l 

Proof. Given the assumptions on d{x), it is clear that (j3.ip holds on all interior regions, and 
so on all open subsets uj d fl. We thus prove that this inequality holds near the boundary 
F. All of our computations below will hold near the boundary F. 

We use an induction argument. The case s = 1 is of course the classical Hardy inequality. 
Let us now assume that the inequality p.ip holds for a given s > 1, and suppose that 

With a denoting the multi-index (ai, 02, as) and using the notation |a| = ai + 0:2 + 03, we 
set 



oai 00:2 QQ!3 

f^Xl ^X2 ^X3 

and for m £ N, 

|a| — m 

A simple computation shows that for m G N, 
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with 

m 
k=0 

for a constant depending on k and to. From the regularity of u, we see that / € Hq{^1). 
Next, with D = D^, we obtain the identity 



Df = Y, C^D'+^-'^u k\{-l)''d'-'' + J2 CsD'-'^u fc!(-l)'=rf"-'=-i(s - k) 

s 



s-l 



+ Yl C^+^D'-^u (fc + ly.i-lfd'^^-^ 

k=0 

= D'+^u s\{-iyd' . (3.3) 
Since / G Hq{U), we deduce from (|3.3I) that for any x £ il with 0:3 G (0, 

f{xi,X2,X3)^{-iysl / D''+'^u{xi,X2,yz) vldyz, 



which by substitution in p.2p yields the identity 



D%-){xi,x2,y3) = j+i ■ 

a x^ 

Hence, a simple upper-bound provides the inequality 

ins^Uw x| ^ , Mxs) Jo^ \D''+^u{xi,x2,y3)\ dys 

where ipi is the piecewise afhne function equal to 1 on [0, and to on [|, 1]. Next, for 
any x £ with 2:3 £ [^,1), we obtain in a similar fashion that 

m ^ , V'2(a;3)X,MD''+iM(a;i,X2,y3)| d2/3 
1^ ^ ^' d(x) ' 

where tp2 is the piecewise afhne function equal to on [0, j] and to 1 on [^,1]. It follows 
that for any x £ 

,| , tpiixs) Jq' \D''+^u{xi,X2,y3)\ dy3+'ilj2{x3) \D''+^u{xi,X2,y3)\ dy3 

1^ (rf)^")' ^ dix) • 

(3.4) 

Setting 

g = tpi{x3)f \D''+'^u{xi,X2,y3)\ dy + tp2{x3) f \D''+'^u{xi,x2,y3)\ dys , 

Jo Jx3 

we notice that g £ _ffQ(il), with 

Ml < C\\D^+'u\\o. 
Therefore, by the classical Hardy inequality, we infer from (I3.4p that 

\\D%^)\l<C\\g\U<C\\D^+'u\\o. (3.5) 
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Since we assumed in our induction process that our generalized Hardy inequality is true at 
order s, we then have that 

which, together with p.Sp . implies that 

II^IL<chiU-fi, 

and thus establishes the property at order s + 1, and concludes the proof. □ 

3.1. K- independent elliptic estimates. In order to obtain estimates for solutions of our 
approximate K-problem (|7.2p defined below in Section [71 which are independent of the 
regularization parameter k , we will need the following Lemma, whose proof can be found 
in Lemma 1, Section 6 of |9]: 

Lemma 3.2. Let k > and g E L°°{0,T; H''{n))) be given, and let f € H^{0,T; H''(n)) 
be such that 

f + Kft=g tn{0,T)xI. 

Then, 

ll/llL°°(o,T;H=(n)) <C max{||/(0)||s, ||.9||L°°(o,r;ff»(n))}- 

In practice, / will usually denote L{V), where L is some nonlinear (possibly degenerate) 
elliptic- type operator and V is some combination of space and time derivatives of ri{t). 

3.2. The embedding of a vifeighted Sobolev space. The derivative loss inherent to this 
degenerate problem is a consequence of the weighted embedding we now describe. 

Using d to denote the distance function to the boundary F, and letting p = 1 or 2, the 
weighted Sobolev space iJ^p(f2), with norm given by J^d{xy {\F{x)\'^ + \DF{x)\-^) dx for 
any F e if^p(il), satisfies the following embedding: 

therefore, there is a constant C > depending only on il and p, such that 

In 

See, for example. Section 8.8 in Kufner [19]. 



F\\U/2<C d{xr{\Fix)\' + \DF{x)\')dx. (3.6) 



4. The Lagrangian vorticity and divergence 

We use the permutation symbol (j2.1l) to write the basic identity regarding the ith com- 
ponent of the curl of a vector field u: 

(curlu)i = eijku'^,j . 

The chain rule shows that 

(curlM(7y)), = e,jkAp'',s , 
Using our definition (|2.3p of the Lagrangian curl operator curl^, we write 

[curl^w]i := e^jkA'jv'',s ■ (4.1) 
Taking the Lagrangian curl of (|1.10p yields the Lagrangian vorticity equation 

EkjiA'^^v^s = , or curl^ w* = . (4.2) 
Similarly, the chain-rule shows that divu(?7) — A^v^^j, and according to (|2.2p . 

div^w = . (4.3) 
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5. Properties of the determinant J, cofactor matrix a, unit normal n, and a 

POLYNOMIAL-TYPE INEQUALITY 

5.1. Differentiating the Jacobian determinant. The following identities will be useful 
to us: 

dJ = a^d-^—^ (horizontal differentiation ) , (5-1) 

dfJ = a^— — (time differentiation using v — rjt) . (5.2) 
ox'' 

5.2. Differentiating the cofactor matrix. Using (|5.ip and (|5.2I) and the fact that a = 
J A, we find that 

da^ = d-^J^^[ala^ - ala^] (horizontal differentiation) , (5.3) 

'^fO'i ^ -Q-;J^^[K'^i - a^a'^l:] (time differentiation using t; = ryt) . (5.4) 

5.3. The Piola identity. It is a fact that the columns of every cofactor matrix are divergence- 
free and satisfy 

a-,fe = 0. (5.5) 
The identity (|5.5p will play a vital role in our energy estimates. (Note that we use the 
notation cofactor for what is commonly termed the adjugate matrix, or the transpose of the 
cofactor.) 

5.4. A geometric identity involving the curl operator. 
Lemma 5.1. 

dtat, al = [curl curl ^^l'^ + (j-'Ka\ ~ a^^^la^ - ~ 

The structure of the right-hand side will be very important to us: the curl structure of 
the first term will be crucially used in order to construct solutions; the second term can be 
made small by virtue of the fact that J~^\a%a!l — ala^a\ — — can be made 

small for short time; the third term is lower-order with respect to the derivative count on v 
and can be made small using the fundamental theorem of calculus. 

Proof of Lemma \5.1i Using the identity ()5.4p , we see that 

Adding and subtracting [S^Sf — 6f6^]Sl, and using the identity curl curl = Z3div— A yields 
the result. □ 

5.5. Geometric identities for the surface ry(i)(r). The vectors ?7,q for a — 1,2 span 
the tangent plane to the surface r{t) — ri{t)(r) in M"^, and 

Ti :— , T2 := -. r , and n 



\V,l\ \V,2\ |77,iX?7,2| 

are the unit tangent and normal vectors, respectively, to T. 

Let gai3 — r],a ■ri,fj denote the induced metric on the surface F; then detg = \ri,i xr],2 P 
so that 

y^n := 77,1 xij,2 , 
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where we will use the notation ^ to mean ^/deFg. 
By definition of the cofactor matrix, the row vector 



?7^l'7^2-^y^l?7^2 



, and ^ : 



It follows that 



(5.6) 



(5.7) 



5.6. A polynomial-type inequality. For a constant A/o > 0, suppose that f{t) > 0, 
i I— >■ f{t) is continuous, and 

fit)<Mo + CtPifit)), (5.8) 

where P denotes a polynomial function, and C is a generic constant. Then for t taken 
sufficiently small, we have the bound 

fit) < 2Mo . 

This type of inequality, which we introduced in f9l, can be viewed as a generalization of 
standard nonlinear Gronwall inequalities. We will make use of this inequality often in our 
subsequent analysis, and also in the form wherein ^/t replaces t. 

6. Trace estimates and the Hodge decomposition elliptic estimates 

The normal trace theorem provides the existence of the normal trace w • of a velocity 
field w € L^(0) with divw e L^{^1) (see, for example, [37] )• For our purposes, the following 
form is most useful: if dw G L'^( fl) with divw € L^(ri), then dw ■ N exists in i/~°-^(r) and 



\\dw ■ iV||^_o.5(r) < C \\dw\\l,^,,^ + \\dWw\\l2^n) 



(6.1) 



for some constant C independent of w. In addition to the normal trace theorem, we have 
the following 

Lemma 6.1. Let dw g L'^{fl) so that curlw e L'^{fl), and let Ti, T2 denote the unit tangent 
vectors on T, so that any vector field u on T can be uniquely written as u^Ta,. Then 



\9w ■ Ta\\jj-o.5(r) ^ ^ 



L2(o) 



|curlw|'^ 



a = 1,2 



for some constant C independent of w. 

See [1] for the proof. Combining (|6.ip and 



\\^w\\H-o■^T) < C ||9w||i2(o) + ||divw||L2(o) + 



L2(0) 



(6.2) 



(6.3) 



for some constant C independent of w. 

The construction of our higher-order energy function is based on the following Hodge- type 
elliptic estimate: 

Proposition 6.2. For an H'' domain n, r > 3, if F e L'^{n;R^) withcmlF e i7"-i(f^; M^), 
divF e H''^^{fl), and F ■ N\r G FF^i{r) for 1 < s <r, then there exists a constant > 
depending only on D, such that 



\\F\\s < C niFllo + II curl^^||,_i + II divF||,_i + \dF ■ N\^^ 
\\F\\s < C (||F||o + II curl^^||,_i + II divF||,_i + J2l=i \dF ■ T„ 
where N denotes the outward unit-normal to T, and Ta are tangent vectors for a = 1, 2 



(6.4) 
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These estimates are well-known and follows from the identity —AF = curlcurlF— DdivF; 
a convenient reference is Taylor |36j . 

7. An asymptotically consistent degenerate parabolic k-approximation of 
the compressible euler equations in vacuum 



In order to construct solutions to ()1.9|) , we will add a specific artificial viscosity term to the 
Euler equations that preserves much of the geometric structure of the Euler equations which 
is so important for our estimates, and which produces a degenerate parabolic approximation, 
which we term the approximate K-problem. 

7.1. Smoothing the initial data. For the purpose of constructing solutions, we will 
smooth the initial velocity field uq. We will also smooth the initial density field po while 
preserving the conditions that p(x) > for x d fl, and that po satisfies (|1.5|) near F. 

For 1? > 0, let < e Cjf'(R'^) denote the standard family of mollifiers with spt(£ii5) C 
B{0,'d), and let fji denote a Sobolev extension operator mapping H''{fl) to H''{M.^) for 
s > 0. _ 

We set Uq — ^ Sniuo), so that for > 0, Wq G C°°(r2). The smoothed initial density 
function is defined as the solution of the fourth-order elliptic equation 

A^p^ = g^*fo(A2po) in (7.1a) 

p^ = on F, (7.1b) 

dpt . dpo , . 

~ A^— — on 1 , (7.1c) 

(xi, X2) I— ^ Po(2;i, X2, xs) is 1-periodic . (7. Id) 

is the boundary convolution operator defined in Section 18.6.11 By elliptic regularity, 
Pq G C°°(ri), and by choosing > sufficiently small, we see that Pq{x) > for a; S $7, 
and that the physical vacuum condition ()1.5|) is satisfied near F. This follows from the fact 



that < on F for 1? > taken sufficiently small, which implies that Pq (x) > ior x £ fl 
very close to F. On the other hand, Pq{x) > for a\\ x E oj for any open subset w C f2 by 
taking sufficiently small. 

Until Section 110.51 for notational convenience, we will denote Uq by uq and p^ by pq. 
In Section 110.51 we will show that Theorem 11.11 holds with the optimal regularity stated 
therein. 

7.2. The degenerate parabolic approximation to the compressible Euler equa- 
tions: the K-problem. 

Definition 7.1 (The approximate K-problem). For k > 0, we consider the following se- 
quence of degenerate parabolic approximate K-problems: 

povl + a'l{plJ-^),k +>idt[a';{plJ-^),k] =0 m ^Ix (0,rj , (7.2a) 

(?7,u) = (e,ito) on 17 x {t 0} , (7.2b) 
po = on F . (7.2c) 

Solutions to (|1.9p will be found in the limit as k — 0. 

Note that (rr2a|) can be equivalently written in a form that is essential for the curl 
estimates that we shall present below: 

i-J-f 2Af(poJ-i),fc+2^9t[A'(PoJ"'),fc] =0. (Eli) 
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Remark 7.2. There appear to be few other possible choices for the artificial viscosity term 
given in \7.2a^ . Our choice, AtSt [af (/r?o^^~^),fe ] , preserves the structure of the energy es- 
timates and also, thanks to Lemma \3.Sl the structure of the elliptic-type estimates that we 
use to bound normal derivatives. On the other hand, the addition of this artificial parabolic 
term does not exactly preserve the transport structure of vorticity, but instead produces error 
terms that we can nevertheless control. 

Remark 7.3. Note that we do not require any compatibility conditions on the initial data in 
order to solve the Euler equations 11. 9\) (or in Eulerian form il.l]) ). and the same remains 
true for our approximate n-problem |7.^[ ^. The lack of compatibility conditions stems from 
the degeneracy condition fil.5\) which allows us to solve for rj and v without prescribing any 
boundary conditions on displacements or velocities. 

7.3. Time-differentiated velocity fields at i = 0. Given uq and po, and using the fact 
that ri[x,Q) = x, the quantity vt\t=o for the degenerate paraboUc K-problem is computed 
using (rr2al ): 

vl\t=o = - {2Kdt[A'l{pQj^'^),k] + 2A'l{poJ-'^),k)\^^f^ = (2KpodivMo - 2po),i +2ku'^„ po,k ■ 
Similarly, for all fc G N, 



afe-l 



dtf" 



- {-2^dt[A'y{poJ),k ] - 2A?(po J-'),fc ) 



These formulae make it clear that each dfV\t=o is a function of space-derivatives of uq and 

7.4. Introduction of the X variable and the K-problem as a function of X. We 

consider a heat- type equation which arises by letting a-'^dx act upon equation (j7.2ap . and 
using the Piola identity (|5.5p : 



aia'^-ipldtJ-%k],j^-4aidta'l-iplJ-'),kh~2[ai^^ . (7.3) 

Po J Po 



Since 5t J"^ = -2J-^Jt, we write ((73)) as 



PO 



1 



-n[aldta'y-ipp-%k],,'2[alA'y{poJ-'),k] 
Po 



Definition 7.4 (The X variable). We set 

X = poJ^^Jt = paJ^^a^v'' ,s = poJ~'^ div^ v . 
Using (|7.5p . we see that 



(7.4) 
(7.5) 



alvl,j = Jtt - dtal 



— 
3 — 



Po 



so that we can rewrite (17. 4p as the following nonlinear heat-type equation for X: 



J^Xt 
Po 



- 2k 



aia'l — ipoX),k],j ^ -K[aldta'l — iplJ ^),k],j -2[aiA'l{poJ ^),k],j 
Po J Po 



It follows from ([73)1 that 



div,, V 



-SJ-'iJt)^+dta^v\, 

__{XJ^^ 
Po 



(7.6) 
(7.7) 
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SO that time-differentiating ()7.7p . we see that 



diY^ vt = ^^^^ ~ dtAiv\, . (7.8) 
Po 

7.5. The nonlinear Lagrangian vorticity equation. The analogue of (|4.ip for our ap- 
proximate K-problem takes the form, with / = poJ^^, 

curl^ vt 2Ke.jiv'',s [/,, A^] AJ" 

^2Ke.j,v\sAt[p,rj{r])], (7.9) 

where / — p{r]). 

We now explain how the formula (|7.9p is obtained. We have that the kth component of 
the Lagrangian curl is 

[curl^ vtf = -2Kekji[dtA\f,i +A\dtf,i],r A] 

= -2nek,AdtA\f,i],rA], 

where we have used the Lagrangian version of the fact that the curl operator annihilates 
the gradient operator; namely ekji{A\F,r ),j = for all differentiable F. 

It is now convenient to switch back to Eulerian variables. We expand dtA\^ and write 

dtA\f,i = -v\s Atf,i Al = p,r ] o r, . 

Now we can compute the standard curl operator of this quantity to find that 

^kji{^ n Pit ] 5j — ^kji'^ ^ij P:r ~^^kji^ :i P:rj 
— ^kji"^ n Pirj • 

Reverting back to Lagrangian variables yields the identity (|7.9p . 

7.6. A boundary identity for the approximate K-problem. For the purposes of con- 
structing solutions to (j7.2l) we will need the formula for the normal (or vertical) component 
of Vt on F: 

= -2J-^alpn,:i~2Kdt[J^^al]p^,^ 

= -2J-^4po,3 ~2Kj-^dtalpo,3 -2KdtJ-^ a^^o.a , (7.10) 

^3 = (^7)1 x^y^a) • 63, (7.11) 

dtal (u,i x?7,2 +?7,i xu,2 ) • 63 . (7.12) 

We note for later use that linearizing (j7.12p about rj — e produces divr v as the linearized 
analogue of dta^. 

8. Solving the parabolic k-problem (17.21) by a fixed-point method 



where 



8.1. Functional framework for the fixed-point scheme and some notational con- 
ventions. For T > 0, wc shall denote by Xt and Yt the following Hilbert spaces: 

Xt = [v e L2(0, T; H\n})\ d^v e L^{0, T; H^-'^ifl)) , a = 1, 2, s} , 
YT = {ye L'{0, T; H'im d^y € ^^(o, T; ij3--(r!)) , a = 1, 2, s} , 
Zt = e Xt I poDv e 1^t| , 
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endowed with their natural Hilbert norms: 

3 3 

"^°"''lL2(o,T;H3-»(r!)) ' 



a=0 a=0 



s.nA\\v\\l^ = \\v\\],^ + \\p,DvfY^. (8.1) 

For M > 0, we define the foUowing closed, bounded, convex subset of Xt'- 

Ct{M)^{v^Zt : \\v\\l^<M}, (8.2) 

where we define the polynomial function A/q of norms of the initial data as follows: 

A/'o = i'(||wo||ioo,||po||ioo). (8.3) 

Since we have smoothed the initial data uq and po, we can use the artificially high H^^^{Vl)- 
norm in Aq. Later, in Section Il0.5( we produce the optimal regularity for this initial data. 
Henceforth, we assume that T > is given such that independently of the choice of 

V e Ct(m), 

ri{x,t)—x+ I v{x,s)ds 
Jo 

is injective for t E [0, T], and that 

13 — 

- < J{x, t) < - for G [0, T] and X e f7 . 

This can be achieved by taking T > sufficiently small: with e{x) = x, notice that 

\\Ji;t) - l|!Loo(n) < C\\Ji;t) - 1||2 = II f al{-,s)v\s {■,s)dsh < CVtM . 

Jo 

In the same fashion, we can take T > small enough to ensure that on [0, T] and for 
some A > 0, 

2X\^\^ <alix,t)a'^ixMj^k y^eR^xefl. (8.4) 

The space Zt will be appropriate for our fixed-point methodology to prove existence of a 
solution to our degenerate parabolic k problem (j7.2p . 

Theorem 8.1 (Solutions to the K-problem). Given smooth initial data with pQ satisfying 
Po{x) > for x €z ^ and verifying the physical vacuum condition lll.5\) near T, for > 
sufficiently small, there exists a unique solution v € Zt^ to the degenerate parabolic k- 
problem {7.2^. 



The remainder of Section [8] will be devoted to the proof of Theorem 18.11 

8.2. Implementation of the fixed-point scheme for the K-problem (|7.2p . Given 
V S Ct{M), we define fj{t) = e + v{t')dt' , and set 

A=[Df]]~'^, J = detDfj, anda = J^. 

Next, we set 

B^'^ — dla^ the positive definite, symmetric coefficient matrix. 
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Linearizing ()7.6p . we define X to be the solution of tlie following linear and degenerate 
parabolic problem: 



Po 



2K\B^'' — ipoX),k],j inf]x(0,T«], (8.5a) 

L Po J 

X^O onrx(0,rj, (8.5b) 

{xi,X2) I— > X{xi, X2, X3,t) is 1-periodic , (8.5c) 
X ^Xq:^ Po div uq on fl x {t ^ 0} , (8.5d) 

where the forcing function G is defined as 

G = -4aidta^-iplJ-'),k],j ~2\aiA^ipoJ~'),k],j -iJ-\Jtf + d^v^,j . (8.6) 
Po 

We shall establish the following 

Proposition 8.2. For T > taken sufficiently small, there exists a unique solution to 18. 5]) 
satisfying 

\\X\\x^ < ^fo + TP{\\v\\l^) + P{\\ cuAvr^^) , 

with the norms Xt, Yt, and Zt defined in \8.1\) . and once again P denotes a generic 
polynomial function of its arguments. ( Generic constants are absorbed by the constants in 
our generic polynomial function P.) 

The proof of Proposition 18.21 will be given in Sections I8.5.1I - I575751 

8.3. The definition of the velocity field v. We will define a linear elliptic system of 
equations for v which should be viewed as the linear analogue of equations (|7.8p . (I7.9p . and 

Definition 8.3 (The linear system for the velocity- field v{t)). With v € Ct{M) given, and 
X obtained by solving the linear problem i8.5\) . we are now in a position to define v{t) on 
[0, Tk] by specifying its divergence and curl in fl, as well as the trace of its normal component 
on the boundary T: 

\XJ^] 

div = div?;t — div^^ vt H ^ j (8.7a) 

Po 

curlwt — curlwt — curl^ vt + 2Ke.jiV,l A* S,:^ {f}) + ^ in Q , (8.7b) 

+ 2kpo,3 divr v = 2kpo,3 divr v - 2po,3 J'^^l 

-2Kpo,3J^^dtal~2Kpo,3dldtJ~^ + c{t)N^ onT, (8.7c) 

/ vtdx = -2 f Ai{poJ-'),kdx-2K I dt[Ai{poJ-^),u]dx, (8.7d) 

{xi,X2) '-^ Vt{xi,X2,X3,t) is \ -periodic , (8.7e) 

where the presence of (Iyvtv in \8. 7c^ represents the linearization of dto^ about fj = e, and 
where 

al = e3- {f],i xfj,2), 
i9ta3 = 63 • (w,i x^,2 +77,1 xw,2 ) , (8.8) 
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the function c{t) ( a constant in x) on the right-hand side of {8. 7c[ ) is defined by 

c{t) =\ [ {di^vt - dh,r,vt)dx + \ [ l^^da; - \ j dt^v\, dx + [ J-^dlpo,3 N^dS 
^ Jn ^ Jn Pa ^ Jn Jr 

+ K j J~^dtalpo,zN^dS + K j dtJ~^alpo,3N^dS + K J divr(v - w)po,3 N^dS , 

(8.9) 

and where the vector field S(7y) on the right-hand side of (ig. 7b\ ) is defined on [0,T] x ^l as 
the solution of the ODE 

Vt + 2E{fj) + 2K[E{fi)]t=0, (8.10a) 
S(0) = Dpo . (8.10b) 

The vector field £ on the right-hand side of US. 7b\ ) is then defined on [0,T] x by 

e =2Al^,j+2K[2Al^p„]t, (8.11) 
where tp is solution of the following time- dependent elliptic-type problem for t G [0, T]: 

2[if ]« +'^l^^t[^'^|:,j ]^ = div(cur^ vt - 2Ke.j^v,l Al S,^ [f])) in Vl , (8.12a) 

= on r, (8.12b) 

{xi,X2) ^i^(xi,X2,X3,t) is 1 -periodic , (8.12c) 

^1(^0 = in (8.12d) 

so that we have the compatibility condition for {8. 7b^ 

div(- curl,-^ Vt + 2Ke.j,v,l E,i (77)) + I) = in n x [0,T] . (8.13) 

An integrating factor provides us with a closed- form solution to the ODE (|8.10p , and by 
employing integration-by-parts in the time integral, we find that 

E{fj){t, •) = e'iDpoi-) - / ^vt{t', ■)dt' , 

= e-^Dpo{-) + / ^vit', ■)dt' - —v{t, ■) + ^-^u^{-) . (8.14) 

The formula ()8.14p shows that S(?7) has the same regularity as v. This gain in regularity 
is remarkable and should be viewed as one of the key reasons that permit us to construct 
solutions to (j7.2p using the linearization (18. 7p with a fixed-point argument. 
Similarly, we notice that 

2[^^ V,j (i, •) = / ^ div(cur^ Vt - 2Ke.j,v,: Af E,i {fj)) {f , •) dt' 







* 1 



[v\r A:v^,i A'^ - -V\r A] - 2kvZ A^ %i (i?)] [f , ■) dt' 



V 



Zk ■' k 

t' -t 

e ^ 

-I — - — div curls w (8.15) 
2k 

Since we can rewrite the left-hand side of (|8.15p as 2A^p -f 2[{Al — Sj)'tp,j and with 
G Ct{M), the elliptic problem (|8.15p is well-defined and together with the boundary 



22 



D. COUTAND AND S. SHKOLLER 



condition (|8.12l b) provides the following estimates for any t e [0,T]: 

\\m\U < A/'o + c(T\\vit)\U + Ma) , (8.16a) 

mt)h < A/'o + c(T\\vt{t)h + \\v{t)h + J* Ms) ■ (8.16b) 
Using ((5?TT]) . the estimates lead to 

||£||2 < AAo + c(T\\vt{t)h + \\v{t)h + Ms) , (8.17a) 

11^ Clla <^fo + c{T\\v{t)\U + ^ II^IU) . (8.17b) 

Remark 8.4. The function c{t) is added to the right-hand side of ^8. 7c[ ) to ensure that the 
solvability condition for the elliptic system |<$. 7| ) is satisfied; in particular, the solvability 
condition is obtained from an application of the divergence theorem to equation |<$. 7ap . 

Remark 8.5. Condition (f^. 7d^ is only necessary because of the periodicity of our domain 
in the directions ei and 62. In particular, our elliptic system is defined modulo a constant 
vector, and the addition of c{t)N^ to the right-hand side of |<$. 7c[ ) fixes the constant in 
the vertical direction, while the condition 7d[ ) fixes the two constants in the tangential 
directions. The particular choice for the average of vf , a — 1,2, permits us to close the 
fixed-point argument, and obtain the unique solution of |7.i?p . 

8.4. Construction of solutions and regularity theory for X and its time deriva- 
tives. This section will be devoted to the proof of Proposition 18.21 

8.5. Smoothing v. We will proceed with a two stage process. First, we smooth v and 
obtain strong solutions to the linear equation (|8.5p in the case that the forcing function 
G and the coefficient matrix T}^^ are C°°(r2)-functions. Second, having strong solutions to 
(|8.5p . with bounds that depend on the smoothing parameter of w, we use nonlinear estimates 
to conclude the proof of Proposition 18.21 

Using the notation of Section \TA\ for each t G [0, T^] and for u > Q,we define 

v^'i-.t) ^ Qy* En{v{-,t)) , 

so that for each v > 0, v''{-,t) e C°°{Q). We define G" by replacing A, a, J, and v in (|8?6| 
with A'^ , ay , J'^, and v'^ , respectively. The quantities A'^, , J" are defined just as their 
unsmoothed analogues from the map ff = e -\- /q "S". We also define [B^}'^ = {ay){{a'^)^; 
according to (|8.4I) . we can choose v > Q sufficiently small so that for t e [0,Tk], 

A|CI' < [S''P(:r,i)Oa V^eE3,xel]. (8.18) 

Until Section [8.5.61 we will use B'^ and G" as the coefhcient matrix and forcing function, 
respectively, but for notational convenience we will not explicitly write the superscript v. 

8.5.1. L?'{Q,T; Hq{VI)) regularity for Xttt- We will use the definition of the constant A > 
given in (|8.4p . 

Definition 8.6 (Weak Solutions of (gSl)). A function X e L'^{Q,T; HI{VL)) with ^ e 
H^^{il) is a weak solution of i8.5\} if 
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(i) for all We H^{n), 



n Po 



(8.19) 



PO 

(ii) X{0) = Xo . 

The duality pairing between Hq (ft) and H~^(fl) is denoted by {■,■) , andG £ L'^(0,T; H~^{rt)). 

Recall that if G G H-^{n), then ||G||h-i(j^) = sup{(G', W) | W G H^iil) , = 
1}. Furthermore, there exist functions Go, Gi, G2, G3 in L^(f2) such that (G, W) = GoyV+ 
GiW,i dx, so that llG'||^_i^f^j — inf X]a=o ll^alloi the infimum being taken over all such func- 
tions Gq. 

Lemma 8.7. If G e L^{0, T; H-'^{n)) and ^ e i^(f^), then for T > taken sufficiently 
small so that {8.^^ holds, there exists a unique weak solution to \8.5\) such that for constants 
Gp > and G„a > 0, 
2 

+ sup G 
L2(o.T;ff-i(n)) te[o,T] 

X, ' 



Xt 



< 



m 


2 


\fpo 





GkA 


G 



+ ll"''^IL2(o,T;ifl(n)) 



lL2(0,T;ff-i(O)) 



Proof. Let (e„)„gN denote a Hilbert basis of 7jg(ri), with each e„ being smooth. Such a 
choice of basis is indeed possible as we can take for instance the eigenfunctions of the Laplace 
operator on with vanishing Dirichlet boundary conditions on F and 1-periodic in e\ and 
62. We then define the Galerkin approximation at order n > 1 of (|8.19p as being under the 
form X„ = X;r=o ^"(Oe* such that: Ml e {0, 



Po 



L2(n) 



Rjfc 

Po 



(Go, e/;)i2(f2)- (Gi,— )i2(Q) in [0,r], 

(8.20a) 
(8.20b) 



A^^(0)==(Xo,e,)i2(n). 



Since each ei is in i7'''+^(ri) n for every fc > 1, we have by our high-order Hardy- type 

inequality p.ip that 



Po 



e i7'=(r2) for fc > 1: 



Furthermore, as the ei are linearly independent, so are the and therefore the deter- 

v^o 



therefore, each integral written in (j8.20p is well-defined. 

Furthermore, as tt 
minant of the matrix 

-*'\/Po' N/Po''^'^"^J('j)eN„={l,...,n} 

is nonzero. This implies that our finite-dimensional Galerkin approximation (|8.20p is a well- 
defined first-order differential system of order n + 1, which therefore has a solution on a time 
interval [0, T„], where T„ a priori depends on the rank n of the Galerkin approximation. In 
order to prove that Tn = T , with T independent of n, we notice that since X„ is a linear 
combination of the eg {£ £ {1, ...,n}), we have that on [0,r„], 



J^Xnt 

Po 



Xr 



2k 



B^'' djpoXn) dXn 
Po dxk ' dxj 



(Go, Xn)L^(n) — {Gi 



dX„ 
dxi 
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Since 



n Po 



2k / {pQXn),k Xn,j da; = 2k / B^^X^.k Xn,j dx + 2K / P(),kXnXn,j dx 



and 



2k I po,kXnXn,j dx = -K 

In Po 



PO,jk 

n Po 



BJ''\Xn\^ dx + K 



n po 



P^±pLB^^\X^\^ dx 
n Po 



n I P^B^'^,^\Xn\' dx, 



n Po 



it follows that on [0,r, 
1 d 



J3\^^dx + 2K [ B^^Xn,k Xn,j dx + K 

Jn 

1 
2 



2 dt Po 



Po,k Po,j 
n Po 



B^''=|X„p dx 



\Xn\ 

Po 



-dx + K f ^B^fc|X„|2 dx + K f P^&\^\Xn\^dx 

Jn Po Jn Po 



GoXndx— / GiXmidx. 



Using (|8.4p . we see that 



1 d 



J 



3 1"^" 



2 dt Po 
1 



dx + 2KX / \DXn\'^ dx + kX 



< 



J^|X„p dx 



n Po 



\\-iJnt + Kpo,jkB^' + Kpo,kB^',,\\L^^a) / -\Xn\^ dx + C\\G\\H-iin)\\DX,,\\o. 

n Po 

(8.21) 



1 



Using the Sobolev embedding theorem and the Cauchy- Young inequality, we see that 



1 d 



2 dt Jq Po 



j3\2^dx + K\ I \DXn\^dx + K\ 



^|X„p dx 
n Po 



1 



<C\\-iJ%+Kpo,jkB^'' +Kpo,kBl'',,\\2 / —\Xn\' dx + C^x\\G\\ij-, 

n po 



i(n) 



where the constant Cka depends inversely on kX. Since v G Ct{M), we have that on [0,r], 
'■* 1 - 

\\-iP)t + mhjk + npo,k hdt < GmVI 

for a constant Cm depending on A/, so that Gronwall's inequality shows that Tn — T (with 
T independent of n e N), and with i < J for all v G Ct{M), we see that 



sup C 


x^{t) 


f 

+ kX 


ie[0,T] 


Vpo 


Jo 



X{0) 


f 


Vpo 


"'0 



Setting Cp = p_f 



we see that 



sup G 


Xn{t) 


2 


^Cp \\x,M\l< 


X{0) 


2 

+ C,x 


te[o,T] 


Vpo 





Jo 


^/Po 


Jo 



i(n) 



Thus, there exists a subsequence {X„„} C {^n} which converges weakly to some X in 
L'^{0,T;H^{n)), which satisfies 

1*^(^)11^-1(0) 



sup G 


X{t) 


2 


r\\x{t)\\i< 


X{0) 


+ C,A / 


te[o,T] 


Vpo 





Jo 


y/PO 


"'0 
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Furthermore, it can also be shown using standard arguments that this X verifies the identity 
()8.19|) . It then immediately follows that 



Po 



and that X{Q) = Xq. Uniqueness follows by letting W = X in (I8TT9)) . □ 
Since 11^11^2(0 ^.^2(5^)) < -P(||w||3cy )i it thus follows from Lemma [5771 and (|8.3p that 



Xt 



Po 



+ sup 

L2(0,T;H-i(f2)) te[0,T] 



II^IIl2(0,T;H1(O)) - 



.22) 



In order to build regularity for X , we construct weak solutions for the time-differentiated 
version of (j8.5p . It is convenient to proceed from the first to third time-differentiated prob- 
lems. We begin with the first time-differentiated version of (|8.5I) : 



2k 



Po 



Po 



{poX 



t),k 



ij - Gt - 

Xt^O 
Xt^X^ 



where the initial condition Xi is given as 



Xi = 2kpo 



{PoXo)n 



po 



the additional forcing term C/i is defined by 



Gi = 2k 



B 



1 

Po 



-{poX),k 



^Giinnx (o,r«], 

onTx (0,TJ, 
on 17 X {t = 0} , 

-poG(O) , 
Po 



(8.23a) 

(8.23b) 
(8.23c) 

(8.24) 
(8.25) 



-^0 — Po div 1*0 , and 

(5(0) = —2k curl curl Mo • Dpo — 2Kdivuo Apo + 2nuQ,i po,ij —2Apo — 2(divuo)^ — Uo;i ■ 

According to the estimate (|8.22p . HGf^i ||^2(o T-H-^{n)) — heiice by Lemma 15771 (with Xt, 
Gt + Gi, Xi replacing X, G, Xq, respectively), 



Xf, 



PO 



+ sup 



Xtit) 



Co \\Xt 



lL2(0.T;Hi(n)) 



< C. 



Next, we consider the second time-differentiated version of (18.51) : 
2k 



PO 



{poXtt),k ,3 = Gtt + ^2 in r! X (0, , 

Po J 

Xtt^Q 



Xf 



x^ 



on Fx (0,r,], 
on 17 X {t = 0} ; 



where the initial condition Xi is given as 



X-i = 2kpo 
the forcing function G2 is defined by 



(po^i 



Po 



G2 = dtGi + 2k 



Bt-{poXt),k 
Po 



^+PoGl{0), 

{■P)tXtt 



Po 



(8.26) 

(8.27a) 

(8.27b) 
(8.27c) 

(8.28) 
(8.29) 
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We do not precise ^i(O), but note that its highest-order terms scale hke either D'^uq or 
PqD^uq or D'^po, so that ||-^/po^/i(0)||o < A/q. Using the estimate (|8.26p . we see that Gu + 
||^/2|||2(o T-H-i(n)) - ^- foUows from Lemma l8Jl (with Xtt, Gu + G2, X2 replacing 

X, G, Xq, respectively), 



Xt. 



Po 



+ sup 

L2(0,T;H-i(f2)) t6[0.T] 



Xttit) 



Finally, we consider the third time-differentiated version of (18.5^ : 



2k 



Po 



Po 



{poXttt),k ,j = Gut + g3mnx (0, , 



Xut = 
Xut — X:i 



onTx (0,rj, 
onnx {t = 0}, 



where the initial condition X3 is given as 



X3 = 2k Pa 
the forcing function is defined by 



(poX: 



2)n 



PO 



G3 = dtG2 + 2k 



Bt-ipoXu),k 
Po 



. 1 



^0^2(0), 

{J%Xut 
Po 



(8.30) 

(8.31a) 

(8.31b) 
(8.31c) 

(8.32) 
(8.33) 



Once again, we do not precise t/2(0), but note that its highest-order terms scale like either 
D'^uq or poD^uq or D^po, so that ||v^^2(0)||o < A/q. Using the estimate ()8.30p . we see that 



\Gut + G3\ 



L2(0,T;_ff-i(O)) 



is bounded by a constant C. (Note that this constant C crucially 



depends on > 0.) We see that Lemma [5771 (with Xut, Gttt + ^^3, X-^ replacing X, G, Xq, 
respectively) yields the following estimate: 



Xt; 



Po 



+ sup 

L2(o,T;ff-i(a)) te[o.T] 



Xtttit) 



(8.34) 



8.5.2. L^{0,T; H^{n)) regularity for Xu- We expand the time-derivative in the definition 
of 02 in ((8?29l) and write 

'Bi^^{poX)J 
L Po ' 



4k 



B: 



t' — {poXt),k 



Po 



-2k 



Po 



Po 



(8.35) 



According to the estimates (|8.34l) . together with the Hardy inequality and the smoothness 
of V, 

dt<G; 



[UtWU 


{JftXtu 


2 

+ 


^{JftXu 


2 


{JfttXt 


D 




Po 





Po 





Po 





hence, (|8.27ap and show that 









du 


KPoX),k 


'J 




L Po J 



< c. 



(8.36) 



L2(0,T;L2(n)) 

The bound (|8.36p together with the fundamental theorem of calculus then provides the 
bound on [0,T]: 

||[— (po^),..]„ \\l + \\dt[—{poX),k],j \\1<G. (8.37) 



Po 



Po 
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From this bound, we will infer that < C and that ||Xt||2 < C, and finally that 

/o ll'^ttlll'^^ — ^- We will begin this analysis by estimating horizontal derivatives of DX. 

Definition 8.8 (Horizontal difference quotients). For h > 0, we set 

a^(.) = "(" + ^^^^")-"(") (a ^1,2), 

and = (9^,5^). 

The variational form of the statement that 

II [ (po^),fc],j II is bounded 

takes the following form: almost everywhere on [0,r] and for /(<) bounded in L^(f7), 

2k I &''X,k <l>,jdx + 2K I &''—X dx ^ I f(t>dx Vcj) e H^{n) . (8.38) 
Jn Jn Pq Jn 

We substitute </) = ~d-^d^X into and using the discrete product rule 

= P'^d^q + dlpq, p\x) = p{x + he^) , 

to find that 

2k [ dl^X,k d^X„ dx + 2K [ d^B^^X,k d^X,j dx 

il 12 

+ 2k [ B^^^'Vo,^ d^(-) B'^X,, dx + 2K [ d^[B^''po,k]- d'^X,, dx 
Jn ^Po^ Jn Po 



13 

n 



f B-^'d^X dx . 



We proceed to the analysis of the integrals ia, a — 1,...,5. By the uniform ellipticity 
condition (|8.18p obtained on our time interval [0,T], we see that 

2K\\\d^DX\\l < il . (8.39) 

The integral 12 can be estimated using the L°°(r2)-i^(fi)- L^(i7) Holder's inequality: 

i2<2K\\B'^&'^\\L^^n)\\Xtuk\\o\\daXu,j\\o 
We then see that using the Cauchy- Young inequality for e > we obtain 

i2<C\\X\\l + 4d^X,,\\l. (8.40) 
The integral is requires us to form an exact derivative and integrate by parts. With 

B^(-)=d^{p^')X'' + Po'diX^--d^po^ + ^, (8.41) 
^Po^ Po Po Po 

we write 13 as 

i,^2K f ^^—f^B^^XB'^X,, dx-2K f ^^^^Vo,^ {-B^po)^B^X„ dx . 

Jn PO Jn PO Pq 
V ^ ' ^ / 

13 a is 6 
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Integration by parts with respect to xj shows that 





r 


dx — K 


Po 


Jn 



( 

Po 



Next, we notice that 



i3a^ = K / B^^-^pQ,kPQ 



d'^X 


f 


dx + K 


Po 


Jn 



jk.h j 



d^X 



Po 



dx 



and therefore, according to (jS.lSp . 

isai > -C'/ikII-D^PoIIl^C'm / 





d'^X 




d'^X 


f 


dx + nX / l-DpoP 


In 


Po 


Jo. 


Po 



dx 



> -ChKC\\d'''X\\^ + kX \Dpo\^ 

Jn 

On the other hand, for e > 0, 

i3a^^<4{B''''''P0l},j\\L^ II 

,d^X , 



d^X 



Po 



Po 



dx 



1lo||5''X,,||o 



(8.42) 



<C\\d^XU + e\\^\\l 

PO 

< C\\d''X\\l + eC\\d''DX\\l, 



(8.43) 



where we have used the Hardy and Poincare inequahties for the last inequahty in (I8.43p . 
Furthermore, 

hb < CW''po,i {-d^Po)\\L-in) ll^llo \\d^X„ llo 
Po Po 



d^Po, 
Po 



ll^lliP^^,.llo 



<C||po||4||^||ip^X,,||o 

<C\\X\\i + e\\d^X„\\l. 



Similarly, we have that 



and finally 



(8.44) 

(8.45) 
(8.46) 

Combining the estimates (I8.39l) - (|8.46p . and taking e > sufficiently small, we find that for 
h > small enough 

\\d'^DX\\l<C{\\Xu\\l + \\fro)<C, 
with C independent of h. It thus follows that 

\\dX\\l<C. (8.47) 

Since po is strictly positive on any open interior subdomain of Q, standard regularity 
theory shows that solutions X (and its time derivatives) are smooth in the interior, and 



U<C\\d^[B^''po,k]\\L-in)\\-\\o\\d^X„\\o 

Po 

<c\\x\\l + 4€^,j\\l, 



k<C\\f\\l + e\\d!},X„\\l 
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hence the equation ()8.5p holds in the classical sense in the interior oi Q. It remains to 
estimate |jX,3 |||2(q T-H'^{n)) ^^'^ ^^^^ pm^pose we expand (^^^^{poX),k ^ ,j to find that 

'X\ _ 1 /rB^'= 

-B-' po,j [ — ],a-B po,a[ — 



^,33 +P0,3 (-) ,3 = ^ ( \—ipoX),k 1 „ -5"'^,3a -S^^, ^,3 -{^"X,^ 

Vpn^ -o33 VL Po -I 

(B^ Vcfc )„ - - ^^>0„ (-) ,a -B"=^PO,a (-) ,3 



(8.48) 
Since 

S"W(^),3-S"^^fc--), 

^ Po ^ Po ^ Po^ 

and since ||po,Q /po||l°°(si) is bounded by a constant thanks to our higher-order Hardy-type 
inequality Lemma [3.11 we see from (|8.47p and ()8.30p that 



X 

We introduce the variable Y defined by 



^,33+P0,3( — ),3|lo<C- (8.49) 



Y{t,xi,X2,X3} ^ r"y3, (8.50) 

Jo Po{xi,X2,y3) 

so that Y vanishes at 2:3 = 0, and will allow us to employ the Poincare inequality with this 
variable. It is easy to see that 

X = por,3 . (8.51) 
Thanks to the standard Hardy inequality, we thus have that 

r,3ll^<qi^ll?<c. 

The estimate (|8.47p then shows that 

\\DY\\l<C, 

and hence by Poincare's inequality, 

\\Y\\l<C, (8.52) 

We notice that 

/ X \ 

X,33 +P0,3 { ) '3 = (Poi^,3 ),33 +P0,3 ^,33 

^PO^ 

= P0^,333 +3po,3 ^j33 +P0j33 Y,3 ■ 

SO that 

l|po^,333 +3po,3i",33llo<t^- 

The product rule then implies that 

l!(pon,333||^<C. 

The same L^(0) can easily be established for the lower-order terms poi^, (poi^),3, and 
(po^))33; for instance the identity (|8.51l) shows that 

UPoY,3),3\\1<C. 

By (|8.52p . we see that ||poi^,33 |lo enjoys the same bound. Since 

(P0>'),33 = Poi^,33 +2po,3 Y,3 +P0,33 Y 
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the estimate (|8.52p proves that ||(poi^),33 ||o ^ C*. By definition of the H^{0, l)-norm, we 
then see that 

/ \\poY{xi,X2,-)\\h3(o i)dxidx2 < C . (8.53) 

Now, thanks to our high-order Hardy-type inequality set on the 1-D domain (0, 1), we 
infer from (18.531) that 

\\Y{xi,X2,-)\\ji2^oi)dxidx2<C. (8.54) 
From (I535|) . 

\\PqY,3 {xi,X2,-) + PO,3Y{xi,X2,-)\\H2(^o^i)dxidX2 < C , 

from which it follows, with (|8.5ip . that 

/ \\X{xi,X2,-)\\H2/oi)dxidx2 < C ; 

hence, 

11^,33 |jo<C^- (8.55) 
Combining the inequalities (|8.55p and (|8.47p . we see that 

\\X\\l<C. (8.56) 

The estimate (I8.56P together with (|8.37p then shows that 

II [ {poXt),k],3 llo is bounded. 

By identically repeating for Xt the i/^(il)-regularity estimates that we just detailed for X, 
we obtain that 

\\Xt\\l<C. (8.57) 
The estimates (|8.56p and (|8.57p together with (|8.36p then prove that 

B^^ - 2 

III {PoXtt),k],j \\L^(O.T;L^(n)) - C ■ 

Po ^ . ^ 

Once again we repeat the estimates for Xtt which we just explained for X, this time L^-in- 
time, and we obtain the desired result; namely, 

II- ^tt|lL2(o,T;ff2(n)) < C*- (8.58) 
It follows that (|8.23p holds almost everywhere. 

8.5.3. L^{0,T-H^{n)) regularity for Xf Using ^2B, we write ([OSal) as 

-2Kdt[—ipoX),kU = Gt- . 
Pq Po 

The estimate (|8.58p together with our higher-order Hardy inequality shows that 

\\dt[ {poX),k],j Wh^o.T.mm <C, (8.59) 

Po \ ■ \ 1/ 

and hence by the fundamental theorem of calculus, 

III {poX),k]„\\l<C. 

Po 
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We employ the identity ()8.4ip to find that 



9" 



PO 



-iPoX) 



PO 



{pod''X),u 



3 

X 
Po 



d''B^''{X,k+Po,k-) 
PO ■ 

^ Po,k jr 

Po Po 



Since the last three term on the right-hand side are bounded in L^(ri) thanks to our higher- 
order Hardy inequality and (|8.56p . we see that 



[ ipod'^X),k 

Po 



<C. 



Now repeating the argument which lead to (|8.56p with d^X replacing X, we find that 



15X11? <C. 



(8.60) 



By differentiating the relation ()8.48p with respect to X3 and using the estimate (|8.60p . 
we see that 



||-''^j333 +P0,3 



PO 



ro<c. 



Now by using the variable Y defined by (|8.50p . we can repeat our argument to find that 
Il-?j333 llo — ^ s-iid hence that 

\\X\\l<C. (8.61) 

From (|8.59[) . we then easily infer that 



jyjk 

II [ ipoXtUli\\ldt<C, 

Po 

so that the argument just given allows us to conclude that 

ll"^tlli2(0,T;H3(O)) - • 



(8.62) 



8.5.4. L^(0, T; i/^(f7)) regularity for X . Bv repeating the argument of Section l8.5.3l we find 
that 

|2 



l^llL2(o,T;ff-4(n)) ^ . 



.63) 



We have thus established existence and regularity of our solution X; however, the bounds 
and time-interval of existence depend on > 0. We next turn to better Sobolev-type 
estimates to establish bounds for X and its time-derivatives which are independent of and 
are useful for our fixed-point scheme. 

8.5.5. Estimates for WXW^^ independent of v . 

Step 1. We begin this section by getting i/- independent energy estimates for the third 
time-differentiated problem (|8.3ip . 

Lemma 8.9. For T > taken sufficiently small and S > 0, 



Xf, 



PO 



+ sup 

L-^(o,T-H-^(a)) te[o,T] 



Xttt{t) 



ll"'^"*llL2(0,T;ffl(n)) 



<Afo+ TP(||X||2 ) + CS\\X\\x^ +TPi\\v\\l) + P{\\ curl5||2 ) . (8.64) 
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Proof. We write the forcing function Gut + Ga as 



Gut + Gs = Gut + dtG2 + 2k 
Ti 



B: 



1 

PO 



iPoXtt),k 



{P)tXut 



(8.65) 



T2 



We test (I8.31a[) with Xut- In the identical fashion that we obtained (|8.2ip . we see that 



1 d 



32 dt po 



dx + 2K\ / \DXut\ dx + nX 



\Dpol 
n Pq 



dx 



< \\^iJ% + Kpo,jk B^^ + Kpa,k B^\, \\L^(^n) / —\Xtttr dx + {Gtu + Gs , X, 
^ Jn Po 

Integrating this inequahty from to f e (0, T], we see that 



— sup / —\Xut{t)\^dx + 2KX 
•J^ te[o,T] Jn Po 



Jn 



\DXut?dxdt<No 



1 

n Po 



\Xttt{t)? dx+ / {Gut + Gs , Xt: 



+ T sup II -(J^)t + KPo,ifc5^''' +K/9o,fcS^'^i 
te[o,T] ^ 

By the Soblev embedding theorem 

\\^{J^)t+Kpo,jk&^ +KPO,kB^'',j h'^in) < G\\^{P)t+Kpo,jk&'' +KP0,k&\3 h- 

The highest-order derivative in the term ^{J^)t scales like Dv, while the highest-order deriv- 
ative in B^'^jj scales like D^ij, which means that we have to be able to bound suptgjg j^j ||w(t)||3 

as well as sup(gjo,T] Il^(^)ll4; ^-^d these are clearly bounded by C\/i||u||xT- Therefore, by 
choosing T sufficiently small and invoking the Poincare inequality, we see that 



C sup 



te[o,T] Jn Po 



1 



\Xut{t)\^dx + C^x 



Jn 



DXutl^ dxdt <Afo+ {Gut + ^3 , Xut)dt . 



We proceed to the analysis of the terms in {Gut + Gz , Xut)dt, and we begin with the 
term Ta in (|8.65l) . We have that 



{rz.Xut)dt< sup ||(P)t||L-(n) 

t£[0,T] 



Jn Po 



\Xut? dxdt<CMT^ sup 



te[o,T] Jn Po 



1 



\Xutf dx. 



where we have made use of the Soblev embedding theorem giving the inequality || {J'^)t\\L°°{n) < 
C||('/^)t||2 < ViCM, where Cm depends on M. 
To estimate the term 72 in (|8.65p . notice that 



Xt, 



{T2,Xut)^~2K / Bi''{Xu,k+Po.k—)Xut,jdx<C\\Bt\\2\\Xu\\i\\Xut\\i 

Po 



and thus 



< 5\\xut\\i + cmiwxuwi < swxutwi + c\mii\\xu{om + twxutimi) , 



{T2,Xttt)dt <Afo + S\\X\\x^ +TP{\\v\\l,) + TP{\\Xfx^) . 



It remains to estimate (7i, Xut)', we use the identity (|8.29p defining G2 to expand Ti as 



Ti = Gtu + dtG2 = Gt 



dttGi 



dt(2K 



Bt-{pat),k 
Po 



{J%Xu 
Po 
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The terms 



{dt(2K 



'■J PO 



, Xttt) are estimated in the same way and 



have the same bounds as {T2,Xttt) and {T3,Xttt) above, so we focus on estimating {Gttt + 
duGi, Xttt) ■ To do so, we use the identity (j8.25p defining Qi and write 

Gttt + dttgi = Gttt ' ^ ] \ ^ ({J%Xt^ 

Expanding Si as 



dtt(2K 



B 



Po 



Ott{^ 



Po 



52 



Si = 2k 



B 



1 



ttt — iPoX),k 
Po 



4k 



B: 



it—ipoXt),k 

Po 



'J 



2k 



B 



t — iPoXtt),k 
Po 



Sit 



Si a 



we see that for (5 > 0, 



{Sia,Xttt} = -2k / Bl^^{X,k+Po,k —)Xttt,]dx 



<C\\Bitt\\o{\\X,k\ 



X 
Po 

,X 



2)\\DXttt\\o<C\\B'i^t\\o\\Xh\\DXttt\\o 



Po 



<c\\Bi^t\\i{\\m\\i+tm\i)+mut\\i. 

where we have used the Sobolev embedding theorem for the first inequahty, our higher- 
order Hardy inequahty Lemma 13.11 for the second inequahty, the Cauchy- Young inequahty 
together with the fundamental theorem of calcuhis for the third inequahty. We see that 

{Sia,Xtu)dt <Uo + S\\X\\xr +TP{\\v\\lr) + TP{\\Xfxr) ■ 

The duahty pairing involving Sn, and Si^, can be estimated in the same way to provide 
the estimate 

'''^{Si,Xttt)dt <Uo + S\\X\\xr +TP{\\v\\%^) + TPiWXWl,^) . 
The duality pairing involving ^2 is estimated in the same manner as Ts and Si to yield 
(52,Xttt)dt <AAo+,5||X|l2 )+TP(||X||2 )+CMri sup / —\Xttt\^ dx. 



te[o,T] Jn Po 

It thus remains to estimate the duality pairing {Gttt, Xttt.dt. We write 
Gttt = - dttthj-\Jt f - dtc4 v\j ) +~2dttt[KAKpoJ-^),k],o ndttt[aldtat-{plJ^^),kh 

^ L ^ V ' Po 



Notice that by the Cauchy-Schwarz inequality 



(£1 , Xttt)dt< r WVP^dttthj-^Jt)' - dtalv\A\\o ll^llo 
Jo ^ ' VPo 



dt 



<TPiMzT) + CT sup / —\Xttt\'dx. 



1 



Next, we write 



{C2,Xt 



te[o,T] Jn Po 
dttt[aiAiiPoJ^^),k] Xttt,j dx 
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The higher-order derivatives in dttt[alA'l;{poJ^^),k] scale hke either D{poDvtt) or Dvtt 
so the fundaniental theorem of calculus and the Cauchy- Young inequality once again shows 
that for (5 > 0, 

i-T 

{C2,Xut)dt = TP{\\v\\l^) + S\\X\\xr- 



A good estimate for £3 requires the curl structure of Lemma 15.11 We write the highest- 
order term in 

Kdut[aidta^ — {plJ-^),k],j = Kdttt[2a{dta^ po,k + poa{dta^ ,k 
Po 



2Kdttt[{aidta'y)po,kJ-%j ; (8.66) 
all of the other terms arising from the distribution of dut are lower-order and can be esti- 
mated in the same way as £2- Now, using Lemma [STTl the highest-order term in (|8.66p is 
written as 

dttt[{aldta'^)po,k J^'^],j = [curl curl i;ttt]''po,fc ^"^ 



J~ 



with TZ being lower-order and once again estimated as £2 • Integration by parts with respect 
to the curl operator in the term ti, we see that 



(il,^t: 



curl wttt • D X {DpoJ ^Xttt)dx 



<\\cuY\vutmo (WDpoJ-^L^in) \\DXtu\\Q + \\cuYl{DpoJ-^)\\LHn) \\Xut\\mn)) 

< II curl z;ut(t) II (||i?/5o J^'lb + || curlppoJ"')||i) ||^m||i 

< C|| cmlvtuiml {\\DpoJ-^\\l + II curl(Z?po J^')||?) + S\\Xttt\\l ■ 



It follows that 
For t2, 

(^2, Xttt) = 



{iuXut)dt < P{\\ cmlvfYr) +TPm\U + S\\X\\x, 



Given that ||a(t) -Id||3 = || at{t')dt'\\3 < VtP{\\v\\zj,), we see that 

{t2,Xut)dt<TP{\\v\\%^) + S\\X\\l,^. 

The duality pairing involving ta can be estimated in the same way. 

Summing together the above inequalities and taking T > sufficiently small concludes 
the proof. □ 
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It is easy to see that we have the same estimates for the weak solutions X, Xt, and Xtt 
solving (|83|) . ([8?23l) . and ([8?27| . respectively: 



3 

E 

a=0 



r' \\f>°- y\\ 



+ sup 

L2(o,T;if-i(f2)) te[0,T] 

< AAo + TP(||X||2, J + C<5||X||2,^ + TP(||C|||J + P(|| curlt;|l|.J . (8.67) 

Step 2. Returning to the definition of Q2 in (|8.29p . by using the estimate (|8.67l) together 
with the Hardy inequality, we see that 

\\Gtt + e2||i.(o,T;L2(n)) < A/'o + TP{\\X\\\^) + C5||X||2,^ +TP{\\v\\l^) + P(|| curlT;|||.J . 
Combining this with the estimate (|8.64p . the equation (|8.27ap shows that 



{PoXtt),k 

Po 



L2(0,T;L2(n)) 



Po 



L2(0,T;L2(a)) 

< AAo + rP(||X||2, J + C<5||X||^^ + TPdl^^lllJ + P{\\ curlvf^^) . 

By repeating our argument of Section 18.5.21 but this time using Soblev-type estimates for 
the horizontal-derivative estimates (replacing the difference quotient estimates as we already 
have regularity), we obtain the desired bound: 



IX; 



tt|lL2(o^T;L2(n)) 



<JVo + TP(\\X\ 



icJ + CS\\X\\t 



-TPi\\v\\%^) + Pi\\cnrlvfY^). (8.68) 



Step 3. From the definition of Qi in (|8.25p . we similarly see that 



IGt + Gil 



L2(0,T;Hi(f2)) 



<Mo+TP(\\X\ 



C6\\X\\lc^+TP{ 



Pdlcurlwl 



Following our argument for the regularity of Xtt, we obtain the estimate 



l^tlli2(o,T;//3(n)) 



2 



,^ <M^ + TP{\\X\ 
Step 4. Finally, 

\\G\\\-2{a,T-m(ii)) < A/'o + TP[\\X\\\^^ 



C5\\X\\\^+TP{ 



C5\\X\\\^+TP{ 



'P(|| curlf 



\Yt> ■ 



.69) 



P(|| curlwll?. 



Following the argument of Step 3 and using Sobolev-type estimates for the horizontal- 
derivative bounds (which replace the horizontal difference-quotient estimates), we finally 
conclude that 

\\ni^^,T-,mm < ^^+TP{m\r)+c5\\nic^ (8.70) 

8.5.6. The proof of Proposition fgTl Summing the inequalities (|8J4| . (fOS]) . (HJH), and 
(|8.70p . we obtain the estimate 

ll^ll^, < AA„ + TP(||X||^ J + C6\\X\\\^ +TP{\\vf^^) + P(|| curli;||^J . 
Choosing 5 > and P > sufficiently small (and readjusting the constants), we see that 

\\X\\\^ < A/'o + TPm\%^) + ^(11 curli;||2.J . 
As the right-hand side does not depend on i/ > 0, we can pass to the limit as ^ ^ in (|8.5p . 
This completes the proof of Proposition 18.21 

Remark 8.10. Suppose that vi and V2 are both elements of Ct{M). For a = 1,2, let Xa 
denote the solution of i8. 5\) with coefficient matrix Ba and forcing function Ga formed from 
Va rather than v. Our proof of Proposition \8.2\ then shows that 

\\Xi - X^fx^ < TP{\\vi - V2\\l^) + P{\\ cmlvi - curlvaHl^^) . (8.71) 

We will make use of this inequality in our iteration scheme below. 
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8.6. Existence of the fixed-point and the proof of Theorem 18.111 

8.6.1. The boundary convolution operator Ae on T. For e > 0, let < € C^(R^) with 
spt(/9e) C -8(0, e) denote a standard family of mollifiers on . With Xh = {xi,X2)^ we 
define the operation of convolution on the boundary as follows: 

KJ{xh) = / p^{xh - yh)f{yh)dxh for / G Ll^{M.^) . 



By standard properties of convolution, there exists a constant C which is independent of 
e, such that for s > 0, 

\K,Fl<C\Fl yPGH^iV). 

Furthermore, 

e\dA,F\o < C\F\o yFeL^{n). (8.72) 

8.6.2. Construction of solutions to |i$.7| j via intermediate e-regularization. We shall estab- 
lish the existence of a solution v to (|8.7p by first considering for any e > the e-regularized 
system, where the higher-in-space order term in (j8.7cp is smoothed via two boundary con- 
volution operators on F: 

div — div Vt — div,^ vt H dtA^v^^j in , (8.73a) 

po 

curluj — curlui — cur^ vt + 2Ke.jiV,l Af E,i {f]) + € in il , (8.73b) 
(u')? -t- 2kpo,3 divr v" = 2kpo,3 divr v - 2po,3 A,[J-'^al] 

-2Kpo,3Ae[J~^dtal]-2KpQ,3A,[aldtJ-^]+c,{t)N'^ on F , (8.73c) 

iv')tdx = -2 f Ai{pQj-^),kdx-2K f dt[Al{poJ~^),k]dx, (8.73d) 

{xi,X2) I— >■ Vt{xi,X2,X3,t) is 1-periodic , (8.73e) 

where the vector 2(77) is defined in (|8.10p and the function Ce{t) (a constant in x) on the 
right-hand side of (|8.73c[) is defined by 

cS) f (div Vt - div,- vt)dx + 1 f I^^dx - ^ / dtAiv\j dx + [ A,[J~^al]po,3 N^dS 
2 Jn 2 Po 2 Jy 

+ kJ^ A,[J-^dtal]po,3 N^dS + kJ^ A,[dtni]po,3 N^dS 

+ K [ divr{A^y - A^v)po,3N^dS . (8.74) 



We now outline the steps remaining in this section. We shall first prove, by a fixed-point 
approach, that for a small time > depending a priori on e, we have the existence of 
a solution to this problem. We shall then prove, via e-independent energy estimates on 
the solutions of (I8.73p . that = T, with T independent of e, and that the sequence v'^ 
converges in an appropriate space to a solution v of (|8.7p . which also satisfies the same 
energy estimates. These estimates will allow us to conclude the existence of a fixed-point 

v = V. 

Step 1: Solutions to (|8.73p via the contraction mapping principle. For 

w e ^ {w e L^{0,T;H^{n)) : wt e L\0,T: H^{n)), wtt e L^{Q,T;H^{n)), 

wttt <^ L^{0,T;H^{n)), (xi,a;2) is 1-periodic}, (8.75) 
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with norm 

II^^IIa'I. = \\'^\\h(o,T;m{a)) + ll^«lli2(o,r;ff3(n)) + ll''"ttllL2(o,T;ff2(n)) + ll^i'ttt|lL2(o,T;/i-i(n)) > 

we set ^{w) = Uo + dt^{w) where dt^{w) is defined by the elhptic system which specifies 
the divergence, curl, and normal trace of the vector field dt^{'w): 

div = div vt - div,-, m + - dtAjv,] in n , (8.76a) 

Po 

curl9t<I>(u)) = curl9tu - curl^ dtv + 2Ke.jiV,l E,l (??) + £ in rj , (8.76b) 
dt^{w) ■ 63 = -2k/9o,3 K divr w + 2k/9o,3 divr v 

-2Kpo,3A4aldtJ^^]+c{w)N^ on T , (8.76c) 

f Ai{poJ-^),kdx-2K f dt[AiipoJ-'),k]dx, (8.76d) 

Vt e [0,T], 9t$(w)(i) is 1-periodic in the directions ei and 62 . (8.76e) 

The function c{w){t) in ()8.76cp is defined by 

If If IXJ'^] If 

[c{w)]{t) ^- / (divut - div^ut)dx+ - / -dx - - / dfAjv^jdx 

^ Jn ^ Jn Po ^ Jn 

+ A,[J-^4]po,3 N'^dS + « ^ A,[J-^dt4]po,3 N'^dS 

+ kJ A,[dtJ-^dl]po,3 N^dS + kJ divr(A2i(; - A,v)po,3 N'^dS , (8.77) 

and is introduced so that the elliptic system (I8.76P satisfies all of the solvability conditions. 
Thus, due to the definition (|8.77p , the problem (18.761) defining dt^{w) is perfectly well- 
posed. Applying Proposition [521 to (|8.76p and its first, second, and third time-differentiated 
versions, we find that 

||at$(,i;)-at$(^S)|U3^ <C(Af,e)||u;--7i|U3^-HCAfr,||«;-^|U3^, (8.78) 

the e dependence in the constant C{M,e) coming from repeated use of (|8.72l) . Note that 
the lack of w on the right-hand sides of (|8.76ap and (j8.76bp implies that both the divergence 
and curl of dt^{w) — dt^{w) vanish, and that on F, 

[(9f$(w) - dt^{w)] ■ 63 = 2kpo,3 Al divr(w - w) + [c(w) - c(w)]N^ . 

It follows from (|8J8)) that 

mw) - Hw)\\x.^ < T,C{M,e)\\w ~ w\\x3^ , 

and therefore the mapping $ : is a contraction if is taken sufficiently small, 

leading to the existence and uniqueness of a fixed-point v"^ = $(w'^), which is therefore a 
solution of (|8J3)) on [0,T,]. 

Step 2: e-independent energy estimates for u*^. Having obtained a unique solution 
to ()8.73p . we now proceed with e— independent estimates on this system. We integrate the 
divergence (|8.73ap and curl (|8.73b[) relations in time, and we now view the PDE for the 
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normal trace (|8.73cp as a parabolic equation for v'^ on F: 



div w*^ — div V — divfj v + 



n 



Po 



curlw^ = curluo + curlu — cur^ v + 2k £-jivZ iv) 



+ / {e.j,v,l dtA'^ + £) in n, 
Jo 

{v')^ + 2kpo,3 divr = +2kpo,3 Ae divr v - 2po,3 AdJ^^al] 



v'{0)=uo, 



.79a) 



(8.79b) 



!.79c) 
.79d) 



[ {v')°'dx = f u^dx - 2 / / 
Jn Jn Jo Ja 



Al{poJ-^),kdx-2K I [A'^;^ipoJ~'),k]dx , 

(8.79e) 
(8.79f) 



yt E [0,T], v'^{t) is 1-periodic in the directions ei and 62 , 

where ?(t) is defined in ()8.74p . 

We will establish the existence of a fixed-point in the metric space Ct„ (M) defined in 
(|8.2p . but to do so we will make use of the space (depending on e) 



^ {w e L^iO,T:H^'^{n))nL'^iO,T;H^{n)) : wt £ L^iO,T; H^{n)), 

wtt e L\0,T;H\n)), wttt e L\0,T;H\n)), A,w G L^{0,T; H*{n)), 
w(0) = Uq} , 

with norm 

II^IIa-I = l|A£w||^2(o,T;H4(a)) + ll'"'t|li2(o,T;//3(n)) + ll'^t«lli2(o,T;H2(n)) + \\wttt\\'l2(^o^T;m{n)) ■ 

Since v S Ct^{M), equations (|8.79ap and (j8.79bp show that both divw' and curlw^ are in 
L2(0, T,; H^{nj); additionally, from (jSjgcl) and (jSTSl) . we see that (i;")^ is in L°°(0, Te; i/^-^(r)), 
and hence according to Proposition 16.21 v'^ G L^(0, T^; i/^(ri)), with a bound that a priori 
depends on e. We next show that, in fact, we can control A^v'^ in CxiM) independently of 
e, on a time interval [0,T] with T > independent of e. 

We proceed by letting d'^ act on each side of (|8.79cp . multiplying this equation by 



P0,3 



d^{v'^)^, and then integrating over F. This yields the following identity: 



-— / \d%vy\'dS + 2n / d-'AidiYrv'd^v'yN^dS 

2dtJrPo,3 Jr 

= [ Gd'^iv'fN^dS - I d'^[po,zKF]d'^{v'f—dS 
Jt Jr Po,3 



where 



G = -2k 



d^po,3 K divr + 39>o,3 dK^ divr + 39^0,3 d^A^ divr w 



F = 2k divr w - 2 J'^Og - 2Kj~'^dtdl - 2KdldtJ'''^ 



(8.81) 
(8.82) 
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Since G contains lower-order terms, we immediately see that for any t £ [OjTg]: 

Gd'^iv'fN^ dS <C\d'^{v')\l . (8.83) 



We then write 

[po,3 A,F] = po,3 d^^,F + d^po,3 A.F + 39^0,3 dA,F + 3dpo,3 d'^KF , (8.84) 

and notice that since the last three term on the right-hand side are lower-order, we easily 

obtain the estimate 

/■_ __ ___Ar3 

I / [(93po,3 A,^^ + 39Vo,3 + 3apo,3 d^A,F)d\v'f dS\ < ClB^'iv^flolB^Dvlo . 

Jr Po,3 

(8.85) 

Next, in order to estimate the highest-order term /p d'^AeFd'^{v^)^N^dS, we notice that 
by the standard properties of the boundary convolution operator Ag, we have that 

J d^A.Fd^iv'fN^dS = J F A,{v')^ dS 

= -2Ky" d^iJ'^dtdl)d^A,{v'fN^dS-2K J d^{dldtJ-'^) d^A,{v'fN^dS 



Jl J2 

+ 2k J d^divrvd^A,{v'')^N^dS-2 J d^[J-^dl] d^A,{v' fN^dS . (8.86) 



1/3 

Jl = -2k / d'\J-''[va xfj,2+fia xv,2?)d'-'A,{v'fN'''dS 



In order to estimate the integral J7i, we recall the formula for dta'^ given in (|8.8I) . and write 



r 

2k I [d^v,i - fj,2iO))Yd'^A,{v'fN'^dS -2k / d'^v\i B'^Aeiv'fN^dS 



J la ^1 fa 



-2k J [(^-^,1(0)) xd'^v,2]'^d'^A,{v'fN'^dS-2K J d'^v\2 B"^ A,{v'f N^^ dS +ni , 

(8.87) 

where TZx is a lower-order integral over F that contains all of the remaining terms from the 
action of B^ ^ so that there are at most three space derivatives on v on F. The trace theorem 
combined with the Cauchy-Schwarz inequality easily show that \R.\\ < C|(w'^)'^|3||ij||4. The 
next crucial observation is that 

Jib + Jid + J3 = 0- (8.88) 
We use the fundamental theorem of calculus, 

51 W - (0) = d,^ and ^{t) - n.i (0) = [dt^JI^ 
to estimate the integrals Jia and Ji^ so that (|8.87p and (|8.88p show that 
\Ji + J3I < Ct\\A,(v'fUv\U + C\{v')%\\v\U . 
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Next, we write the integral J2 as 



where 



7^2-y" d^DT]d^A,{v')^N^dS 



where the symbol ~ is used here to mean that TZ2 is comprised of integrands which have 
the derivative count of the integrand d"^ Dfj A^{v'^)^ . It follows that 

|7^2| < \Dfj\2.5\Mv'fk5 < C||77||4||A,z;l4, 

the last inequality following from the trace theorem. Since f]{t) — e + v, we see that for 
some S > 0, 

|7^2| <Afo + {S + ct^)\\A,vm + ct^Ml , 

Returning to the remaining term in J72, we write 

4k J alJ-^a'^d^v'',s A,{v')^ dS ^ Ak J 8^ dWvB^ A.iv" fN^ dS 



J2^ 

+ 4k / {alJ-^d^.-S'^)d^v\s d^A,{v')^N^dS 



Jib 

We estimate the integral J72a with the Cauchy-Schwarz inequality. The integral J^2b can 
be estimated using the fundamental theorem of calculus: 

\J2b\<Ct\\A,{v'fU\v\U, 

so that we have established the following estimate: 

IJ2I + \Ji + J3I < AAo + TCm + [5 + Ct')\\A,v'\\l + Ct'WvWl 

+ Ct\\K{v'fU\v\U + C\{v'f\:i\\v\U + C\\X\\l + C\\ divi;||3||A,«^||4 , 

where the bound on the integral J^2a in contained in the right-hand side. Finally, the integral 
J74 can be estimated in the same way as 72-2 above, so that with the identities (j8.86p . (|8.85p . 
and (|8.84l) we have shown that 

a3(po,3 A,F)d\v'f dS\ <Ua + TCm + {5 + Ce)\\A,v% + Ct^Ml 

r Po,3 

+ Ct\\A,{v'fU\v\U + C|(«^)3|3||«||4 + C\\X\\l + C\\ divz;||3||A,t;^||4 ■ (8.89) 



We now turn our attention to the second term on the left-hand side of (|8.80l) . which will 
give us as a sign-definite energy term plus a small perturbation. We first notice that by the 
properties of the boundary convolution operator A^, 

d^[Al[v\\+v\l )\d^{v'fN'^dS= [ d^A,{v',l+v'',l)d^A,v' -N dS . (8.90) 



r2 



The divergence theorem applied to the integral X (as our domain fl = T x (0,1)) implies 
that 

I = / B^'A^iv^l +v%l ) B^A,{v')\s dx+ I B^ [Ae(«M3 Ja WA,{v'fdx . 
Jn Jn 
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Now, 



and 



li= - [ \d^A,{v'f,3\^dx+ [ A,d^divv'd^A,v\ldx 



from which it follows that 

I = - [ \d^A,D{v')^fdx+ [ Aj^divv''d^A,{v'f,3dx 
Jn Jn 

+ [ AJ^[cut\v' ■e2]d^A,{v')^,idx'- [ Aj^[cnT\v'-ei]d^A,{v' f,2dx. (8.91) 
Jn Jn 

Now, thanks to (|8.79ap . we have for all t e [0, T^] 

II divi;'|l2 <Ct'\\v\\l + C\\v\\l + C\\^J\l + Ct\\v\\l^oMHHn)) 

<Ct^v\\l + C\\v\\l + Ct\\v\\U,,t,HHn)) + C\\X\\l , (8-92) 

where we have used our higher-order Hardy inequality Lemma |3. II for the second inequality. 
Next, with (I8.79bp . we see that for all t € [0, T,] 

II cuAv'h <Ct\\v\U + C\\vh + C\\uo\U + CVi\\DiEifj))\\ + CVi\\ 

<Ct\\v\\4 + ailuolU + CVt\\v\\L^o,t-HHn)) > (8-93) 

where we have used (|8.17p and the identity (|8.14p . relating S(?7) to v and where we have 
relied crucially on the chain-rule which shows that 

Note that (j8.14p provides us with a bound which is e-independent, but which indeed depends 
on K. 

The action of the boundary convolution operator does not affect these estimates; thus, 
using Proposition 18.21 we see that 

r II dWA^v'Wldt < No + TP{\\v\\%^) + P{\\ curlt;||^J , (8.94) 
Jo 

f ||curlA,«'|j2dt <A/-o+TP(||f;|||J. (8.95) 
Jo 

We integrate the inequality ()8.80p from to i, and we use the estimates (|8.89p . (|8.83p . 
(IHUl), (|5:g5|) together with the fact that ^ > C > on T, to obtain that for any 

te[0,T], 

\d\v^nt)\l + f \\d'DA,{v^f\\l < AAo + Ct\{v^f\l + Ct||i;|||, + Ct f \\A,{v^f\\l 
Jo JO 

+ CVtJ^ Ml + ^J^ \{v')'\l + S\\A,v^\\l + TP{\\v\\%^) + P{\\ curlSlll, J + C|| divf;|||.^ , 

where we have used the Cauchy- Young inequality ab < -^a^ + VtlP' for a, 6 > 0. By taking 
(5 > sufficiently small, and using the relations (|8.94p and (|8.95p . this implies the following 



42 



D. COUTAND AND S. SHKOLLER 



inequality: 

\v%t)\l+ I \\A,v^\\l<Mo + Ct\{v^f\l + Ct\\v\\l+Ct I \\Uv^)Y4 + CVt I \\v\ 



t 

2 

4 



c 



t 

<J\fa + C{t + Vt) sup \v'\l + C(t + Vi)\\v\\l^+Ct I \\Aev'\\i 



te[o,T] 

+ TPi\\v\\%r) + PiW curlz^ll^J + C\\ dWvfY^ . 

(Note that due to the presence of the convolution operators in the definition of (v^)^ on 
r in formula (I8.79cp . it is clear that supjg^t] Iv'^ls is bounded by some finite number, which 
a priori depends on e.) 

Since we are considering a bounded time interval, < i < Cy/t, and we will henceforth 
make use of this fact. The previous estimate then implies that 



sup 

te[o,T] 



\v%t)\l + f WK^v^dt < AAo + VtP{\\v%IA + VtP{\\v\\1^) 
Jo 

+ PiW curl«||2.^) + C\\ diYvWl^ . (8.96) 



Estimates for the time-differentiated quantities are, in fact, very straightforward at this 
stage. By using the expression (|8.79cp . we see that thanks to the estimate (|8.96p . 



Jo Jo 



<Afo + VTPi\\v%.) + VtP{\\v\\1^) + P{\\ curl«||2,J + C\\ diYvWl.^ . 

(8.97) 

Now, just as we estimated the divergence and curl of w'^ in (|8.92p and (|8.93p . we can repeat 
this procedure to estimate the divergence and curl of v^. By using (|8.73ap and (|8.73bp . we 
also have that 

r II dWvtWldt <No + VtP{\\v\\1^) + P{\\ cm\v\\Y^) + C\\ dWv\\l.^ , 
Jo 

II curl^;,'||2di < AAo + VtP{\\v\\%J + P{\\ curl«|||.J + C|| div5|||,^ , 
which in addition to the normal trace estimate (I8.97P provides the estimate: 

T 

e||2 J. ^ «/- , .AFD/IL.e||2 \ , ./^r,i\\r-.\\2 \ , o/||„„ i-||2 n , /^|| j-.-||2 



UUdt <J\fo + VTPiWv^W'^,) + VT P{\\v\\'z^) + P{\\ curlt;||^ J + C|| div«||^^ 

(8.98) 

We proceed in a similar fashion to estimate v^f, by considering the time-differentiated 
version of (j8.73p . and using (I8.96P and (|8.98p . This yields the following inequality: 



/ \\vl,\\ldt <Afo + VtP{\\v%,) + VtP{\\v\\%^) + P{\\ cmlvfYr) + C\\ divv\\ 
Jo 

(8.99) 
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Finally, by using the second time-differentiated version of (|8.73|) and using (|8.96|) . (|8.98|) 
and (|8.99p . we also have that 

r M,,\\ldt <Uo + VTPi\\v%,) + VtP{\\v\\1^) + P{\\ cmlvW^Y^) + C\\ divvWy^ . 
Jo 

(8.100) 

The estimate together with ^Ml, (IH3S1), and (|8.100p provide us with 

\\v'\\%. < AAo + VfPiWvm,) + VtP{\\v\\1^) + P{\\ cm\v\\i.^) + C\\ divz;||^^ , (8.101) 

where the polynomial functions P on the right-hand side are independent of e. 

Thanks to our polynomial estimates in Section we infer from (I8.10ip the existence of 
T > (which is independent of e) such that u*^ G and satisfies the estimate: 

\\v'\\l, < 2Mo + 2Vt Pi\\v\\l^) + P(|| curl5||2,^) + 2C\\ diYvf^^ . (8.102) 

(We will readjust the constant C and the polynomial functions P to absorb the multiplication 
by 2.) 

Step 3: The limit as e ^ and the fixed-point of the map u n- u 

We set e = n g N = {1,2,3,...}. From ()8.102p . there exists a subsequence (still 
denoted by e) and a vector field v&X^,V e L^{0,T; H^{n)) n L°°{0,T; H^{n)) such that 

v*^ ^ V in A"! , (8.103a) 

in A"! , (8.103b) 

A^v" V in L^{0, T; H^in)) . (8.103c) 

where the space is defined in ()8.75p and 

X^ = {weL\0,T;H^{fl)) : wt e L^iO,T; H^{fl)), wu e L''{0,T; H^Q)), 

{xi,X2) I— > w is 1-periodic} . 

Next, we notice that for any ip G C^{n), the space of smooth functions with compact 
support in fi, we have for each i = 1, 2, 3 and t G [0,T], T still depending on k > 0, that 

lim / A,{vy ■ (pdx = \im [ {vy ■ A^(pdx = [ v'lpdx, (8.104) 

where we used the fact that A^ip — > in L^{il). This shows us that v = V, and that 
\\vfx^ <Afo + VfPiMl^) + P(|| cnvlvW^Yr) + C\\ dWy^Y^ . 

The estimates weighted by po i^i the definition of the .ZT^-norm follow immediately from 
multiplication by pq of the equations (|8.7ap and (j8.7bp ; because po vanishes on F and using 
the (unweighted) estimates already obtained, there is no need to consider the parabolic 
equation (|8.7cp . so that 

IkllL < + VtP{\\v\\1^) + P(|l curlz;|l^J + C\\ divf;||^^ . (8.105) 

Moreover, the convergence in (|8.103p and the definition of the sequence of problems (|8.79l) 
easily show us that w is a solution of the problem (|8.7I) : furthermore, we see that we can 
obtain for the system (|8.7I) the same type of energy estimates as in Step 2 above. This 
shows the uniqueness of the solution v of (|8.7p . and hence allows us to define a mapping 
@ : V G Zt ^ V e Zt- 

We next launch an iteration scheme. We choose any w^^^ e Ct{M) and define for n G N, 
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For each n £ N we set r]'^''\x,t) + (x, t')'^^', ^^"^ = ^'"^ = det£ir;("), 

a(n) ^ j(")yi(n)^ is the solution to (IH3]) with a("), and yl^") replacing v, a, 

J, and A, respectively. Similarly, we define S'^"^(?7(")) via equation (|8.14p with w'"^ replacing 
V] we define via equations (|8.11l) and (|8.15p with w*^") replacing v. 
According to (jS.lOSp . for 

< AAo + VTP(||t;(")|||J + P(|| curl«(")||2.J + C\\ div«(")||2.^ . (8.106) 

From (|57a|) . 

divt;(") = div - div„(„-i, + 

so that 



II divz;(")||2.^ < II div«("^i) - div,(„-„ «("-i)|||.^ + I 



Po 

< AAo + VrP(||«("-^)||L) + P(|| curlz;("-i)||2 ) , (g.io?) 



where we have used our higher-order Hardy inequality Lemma 13.11 and Proposition 18.21 for 
the second inequality. Next, we use (j8.7bp and write 

curl?;^"^ = curl Mo + curlw*^""^) - curl^(„-i) u^""^^ 

Jo Jo Jo 

It then follows, using (|8TT4| and ([STH) . that 

||curlM(")|||.^ <AAo + VtP(||i;("-i)|||^). (8.108) 

Combining the estimates (I8.106p . (I8.107p . and (|8.108p . we obtain the inequality 

||^(n+i)||2^^ < AAo + VtP(||i;(")|||J + VrF(||«("-i)|||J + x/TP(||z;("-2) 

This shows that by choosing T > sufficiently small and M >> Ao sufficiently large, the 
convex set Ct{M) is stable under the action of 9. 

In order to see that & has a fixed-point, we simply notice that by proceeding in a similar 
fashion as in Step 2 above (for the e- independent energy estimates), and by using the 
inequality (|8.71l) 

_ < p(||^;(«) _ + Vt P(||i;(«-i) - i;("-')|||^) 

+ Vt p(|| _ ^(„-3) iii^) ^ (8_^Qg) 

where the polynomial function P can be chosen under the form P{z) = X^Jli ^j-^'' ^'^^ some 
integer m > 1. 

Although, the inequality (|8.109l) is not exactly the usual hypothesis of the contraction 
mapping theorem, the identical argument shows that for T = taken sufficiently small, 
the map 8 is a contraction, and possesses a unique fixed-point v satisfying v = Q{v). We 
will next prove that this unique fixed-point v is the unique solution of the K-problem (|7.2p . 
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8.6.3. The fixed-point of the map v ^ v is a solution of the n-problem. In a straightforward 
manner, we deduce from ()8.7|) the following relations for our fixed-point v — v: 



in r2 , 


(8.110a) 


in n , 


(8.110b) 


on r , 


(8.110c) 




(8.110d) 




(8.110e) 



div^v,= ^-^^-dtAlv\, 
Po 

curl^ vt = 2Ke.J^^;,^ Af E,^ (r/) -I- £ 

<dx = -2 [ Al{p^J~^),kdx~2K f dt[Al{poJ-'),k]dx, 
n Jn Jn 

{xi,X2) H> vt(xi,X2,X3,t) is 1-periodic , 
where X is a solution of (|7.6p and where the function c{t) in (|8.110cp is defined by 

~2c{t)= [ {-divvt+div„vt)dx~ [ ^^'^ dx + / dtAjv,)dx-2 [ J^'^alN^ po,3 dS 
Jn Jn po Jn Jr 

- 2k J^J-^dtalN3Po,3dS ~2k J^aldtJ-^N3po,3dS 

divvt + div,^ Vt)dx — / dx + / dtAlv,^jdx 

n Jn Po Jn 

-2 j Dp{ri) ■NdS~2n j [Dp[ri)]t ■ N dS , 

where Dp{r]) = A'^{poJ~^),k, and dS = dxidx2- By using (jS.llOap and the divergence 
theorem, we therefore obtain the identity 

c{t) ^^J^[vt + 2Dp{r^) + 2K[Dp{T])]t]-NdS. (8.111) 

The fixed-point of the map v i-^ v (which we are labeling v as well) also satisfies the equation 

Vt + 2E{T]) + 2K[E{Tj)]t = 0, (8.112a) 
S(0) = Dpa. (8.112b) 

It is thus clear from (|8.11Gcp and (|8.112p that the fixed-point is a solution to the K-problcm 
()7.2p . if we can prove that 

c{t) = and E = Dp, (8.113) 

where we remind the reader that 

p(77) = poJ~' . (8.114) 
This is in fact the case, and we now explain why (|8.113l) holds. 

Step 1. We first apply curl,, to each term of (|8.112ap and compare the resulting equation 
with (|8.110bp . This implies that 

curl^[S(7])] +K[curl„[St(7?)] +u,(77)curl,,[(9,S)(77)]] = -^C, (8.115) 

which thanks to the fact that by definition of curl,,, 

curl,,[e(r;)] = [curie](?7), (8.116) 
for any vector field Q, provides us with 

[curlS](7?) + K[[cnvlEtM + u,{fj)[cur\{d,E)M] = -hi, (8.117) 



46 



D. COUTAND AND S. SHKOLLER 



which shows that 

[curlS](r/) + K[curlS(77)]t = -~D^,{ri) ~ K[D4>,{ri)\^, (8.118) 
where ipe denote the Eulerian version of i/', given by 

V'eO?? = '0- (8.119) 

According to (I8.112bp . S(0) = DpQ\ thus, we have that [curl S] (77) (0) = 0. Furthermore, by 
our definition (|8.12l) . we have Dil)f.{rj)\t=o — in il, which with (I8.118P aUows us to conclude 
that for t e [0,T], 

[curlS + Di^e](?7)(i) =0. 
We may therefore consider the following elliptic problem: 

AV'e = -div(curlS) = in ri{t){n) , 

V'e = on ?7W(r), 

{xi,X2) 'ipeixi,X2,X3,t) is l-pcriodic , 

which shows that V'e = and hence € — 0. Therefore, curlS = in 7y(t, fl) and there exists 
a scalar function Y(t, •) defined on 7/(t, i7) such that 

S = DY. (8.120) 

It remains to establish that DY = Dp. We will first prove that a Neumann-type boundary 
condition plus a small tangential perturbation holds for Y — p; namely, we will show that 
(r-p),3=0onr;(i,r). 

Step 2. We take the scalar product of (I8.112ap with 63 to find that 

+ 2DY{t]) ■ 63 + 2K[DY){T])]t • 63 = 0, 
which, by comparison with (j8.110c|) . yields the following identity on F: 

2 \d{Y - p)ir,) + 2k[D{Y - p)(77)]t] • 63 - -c{t)N^ 

D{~p + Y){tj) + K[Di-p + Y)(r])]t\ ■ N dS N\ (8.121) 



where we have used the expression (|8.11ip for c(t). 
Therefore, since N — (0, 0, N^) on F, this implies: 

[i?(-p+y)(ry) + /t[i5(-p+r)(r;)]t].iVd^= i^[i?(-p+y)(r;) + «[i5(-p+r)(r?)]t].7Vd5, 
and thus that 

-c(t) = J[D{-p + Y)irj) + fi[D{-p + Y){r])]t]-NdS ^0. (8.122) 
The identity (|8.122p together with (|8.12ip then implies that 

D{Y - p){ri) + n[D{Y - p){rj)]t'^ • eg - on F . 

Therefore, 

[{Y - p),3 (77)] {t, •) - e-i{Y - po),3 (0) on F, 
which with (I8.120p and S(0) = DpQ shows that 

(r-p),3=0, onr,{t,T) . (8.123) 
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Step 3. Wc now apply div^ to (|8.112a|) . and compare the resulting equation with ()8.5a|) . 
Using (|8.110ap and the fact that X{0) = podivuo, we have that 

X ^ poJ^^divr,v. (8.124) 

By denoting 

q = p-Y, (8.125) 

this leads us to: 

div^[Dq{r]) + nDqtiv) + KU,{r]){Dq),, {f])] =0 in 17. (8.126) 
This is equivalent in [0, T] x to: 

Aqiij) + KAqt{f]) + KU,{r])Aq„ (r?) + Ajv,) q,u {v) ^ , 

or equivalently, 

Ag(77) + K[Aq{7^)]t + Ajv^j q,u (v) = ■ (8-127) 
Now since po = Y{0), we have that 

Aq{0) =0 infl. (8.128) 
Also, from ()8.123p . we have the perturbed Neumann boundary condition 

(7,3 = 0, onr/(r). (8.129) 
By (|8.110dp . ()8.112p . we obtain that for a = 1, 2, 



q,a (ri)dx + K I [q,a {ri)]tdx = 

Jn 

or equivalently, 

q,a {■n)dx -f ndt / q,a {■n)dx = , 



which together with the initial condition / q,^ {0)dx = implies that 

Jo 



q,c{r])dx^O. (8.130) 

I 

Therefore, by setting / — Aq, we have for all t e [0, T] the system: 

Aq = f in 7j{n) , (8.131a) 

/ J~^q,adx^O, (8.131b) 

q,3 = on ry(r) , (8.131c) 

Dq is 1-periodic in the directions ei and 62 . (8.131d) 

Note that because of the periodicity of u, the domain r/(f2) is such that ry(l,X2,a;3) = 
rjiO, X2,X3) + (1, 0, 0), and r/{xi, 1, X3) = r/{xi, 0, X3) + (0, 1, 0), which explains why condition 
(I8.131dl) holds. 

We now take the vertical derivative ^3 of (|8.131ap . multiply the resulting equation by 5,3 
and integrate by parts in r]{n), using the conditions (|8.131cp and (j8.131dl) . This yields: 

/ \Dq,3\'^dx = - f q,33dx, 

which provides us with the estimate 

\\Dq,3\\oM^)<C\\f\\o,„ia), (8.132) 
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where we are using the notation || • ||o,J7(o) = II • ||L2(r((si))- 
We now write (|8.131ap under the form 

q,ii +q,22 = 9 where g := -9,33 +/ . (8.133) 

It follows from (I8.133P that 

(9,11 +g,22 ) (g,!! +9,22 / \g\'^dx. 

ri{n) Jri{n) 

Integration by parts on the left-hand side of this equation, together with the periodicity of 
Dq and its derivatives, shows that 



q,ai3 q,ap dx + / q,iiq,2n2{t)dS{t) - I q,i2q,2ni{t)dS{t) = I \g\ dx , 

(8.134) 

where the notation dS{t) denotes the naturally induced surface measure on 77(r). Note that 
the presence of these surface integrals is due to the fact that ri{T) is no longer necessarily 
horizontal for t > 0, so that integration by parts in purely horizontal directions produces 
boundary contributions. 
We now notice that 

qaiqan2{t)dS{t) ^ / q,ii q,2 ^-7-rn3{t)dS{t) , 
jj(r) JriiT) "3irj 

where division by n^it) is bounded since for t taken sufficiently small, 1^3 (t)| is very close 
to 1 on ?7(r). 

Next, for i = 1, 2, 3, we smoothly extend ni{t) into ?/(fi), and note the integration-by-parts 
(with respect to X3) identity 

f Ti (i) t Tt (i) 

qaiq,2n2it)dS{t) ^ (7,113 q,2 ^7— + / 9,11 ^3 [9,2 dx . (8.135) 

r,(r) Jriin) n3(t) J,,(f^) n3{t) 

An integration by parts with respect to the variable xi for the first integral on the right- 
hand side of (|8.135p then yields: 



/ 9,11 9,2 "^2 {t)dS{t) = - q,,,d^[q,2^]dx + 9,13 9,2 — TTni{t)dS{t) 



/ 9,11539,2-777 



dx. (8.136) 



Similarly, by integrating by parts first with respect to X3 and then with respect to X2, we 
see that the third integral on the left-hand side of (|8.134p can be written as 



q,i2q,2ni{t)dS{t) = - [ 9,13 52 [9,2 —777] da; + / 9,13 9,2 — 477^2 (i)d'S'(i) 



9,12 93 [9,2 ^]da;, (8.137) 

jj(f2) '^3ltj 
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which shows that the boundary integrals over 77 (F) cancel each other when we substitute 
(|8.137p and (|8.136p into (|8.134p : thus, (|8.134p takes the following form: 




which thanks to the estimate (I8.132P and the relations |»T-a(i)|vi/i'°°(n) ^ Ct, implies that 




g,Q/3 g,a/3 dx < C||/||g_^(jj) + Ct||D2(7||^_^(jj) + Ct\\D'^q\\o,ri(n}\\Dq\\Q^^(n) , (8.138) 



Combining this estimate with (|8.132p . we obtain that 

\\D\\\l,(n) < C|l/llo,,(o) + Ct\\D'q\\l,^n) + Ct\\DW\o^^^n)\\Dq\\oMn) ■ (8.139) 

Now, we notice that the conditions (I8.131cp . (j8.131dp . and (|8.130p yield Poincare inequalities 
for q,a and 5,3, so that 

\\q,3 ||o,,,(f2) < C\\Dq,3 ||o,^(o) , (8.140a) 
< Cp9,a llo.rKO) , a =1,2. (8.140b) 
Therefore, (|8.139p and (I8.140p provide us with 

\\mlnin)<C\\f\\l,^n) + Ct\\Dq\\l^^n^, 
which, by taking T > small enough, yields: 

\\Dq\\lMn) < C\\Aq\\l^^^y (8.141) 

We thus have proved that 

q,u{7i) = F'\t,Aq{7j)), (8.142) 

where i^'*(t, •) denotes a linear and continuous operator from L^(ri) into itself, whose norm 
depends in a smooth manner on v in L^{0, T; H^(n}). 
Therefore, the ODE 

Ag(77) + 4^qiv)]t + Aiv,) F'\t, Aq{rj}) = 0, (8.143) 
with the initial condition (|8.128p allow us to conclude by the Gronwall inequality that on 

Aq{r]) = . (8.144) 

From (|8.144l) and (|8.14ip . we infer that 

Dq = Q in [0,T] x f7, 

which finally proves that DY = Dp, and therefore that S = Dp. Together with (|8.122p . 
this establishes that w is a solution of the K-problem (|7.2p on a time interval [0,rK]. This 
concludes the proof of Theorem 18.11 

By considering more time-derivatives in our analysis, it is easy to see that as long as the 
initial data is smooth, we can construct solutions which are smooth in both space and time. 
We state this as the following 

Theorem 8.11 (Smooth solutions to the k- problem). Given smooth initial data with po 
satisfying po{x) > for x G Q and verifying the physical vacuum condition il.5\) near F, for 
Tk > sufficiently small, there exists a unique smooth solutions to the degenerate parabolic 
K-problem |7. 
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9. K-INDEPENDENT ESTIMATES FOR (17. 2p AND SOLUTIONS TO THE COMPRESSIBLE 

EULER EQUATIONS (|1.9p 



In this section, we obtain estimates for the smooth solutions to (|7.2p . provided by Theorem 
18.111 whose bounds and time interval of existence are independent of the artificial viscosity 
parameter n. This permits us to consider the limit of this sequence of solutions as k — >■ 0. 
We prove that this limit exists, and that it is the unique solution of ()1.9|) . 



Notation. For the remainder of SectionlH rj{t) denotes the solution of the K-problem (|7.2|) 
on the time interval [0,Tk]. In particular, ri(t) is an element of a sequence of solutions 
parameterized by k > 0, but in order to reduce the number of subscripts and superscripts 
that appear, we will not make this sequential dependence explicit. The reader should bear 
in mind that rj is really 

9.1. A continuous-in-time energy function appropriate for the asymptotic pro- 
cess K 0. 



Definition 9.1. W^e sef on [0, T„] 

4 4 



a=0 
4 

E 

a=0 



[II V^po a^a4-'^i?«(s)||2] ds + II curl,, v{t)\\l + IIpo curl,, v{t)\\l . (9.1) 



The function E{t) is the higher-order energy function appropriate for obtaining K-independent 
estimates for the degenerate parabolic approximation {7.0^ . 

According to Theorem 18. Ill solutions to our approximate K-problem (17.21) are smooth, 
and hence T i— >■ supjgjQ j^j E{t) is a continuous function on [0,Tk] to which the polynomial- 
type inequality of Section [5?6l can be applied. 

Definition 9.2. For the remainder of the paper, we will use the constant AIq to be a 
polynomial function of E{0) so that 

Mo = P(^(0)) . (9.2) 

Remark 9.3. Note the presence of n-dependent coefficients in E{t) that indeed arise as a 
necessity for our asymptotic study. The corresponding terms, without the k, would of course 
not be asymptotically controlled. 

Remark 9.4. The 1 is added to the norm to ensure that E{t) > 1, which will sometimes 
be convenient, whereas not necessary. 

Remark 9.5. Of all the terms on the right-hand side of i9. the sum J2t=o[\\Bt"'''l(t)\\1-a 
is the fundamental contribution, providing the basic regularity for the solution. Notice that 
only the even time derivatives of rj{t) appear in this norm. While it is possible to also 
obtain estimates for the odd time derivative of r]{t), we will instead rely on the following 
interpolation estimate: For k > 1, given a vector field r £ L°°([0, T]; _ff'^(n)) with rtt G 
L°°{[0,T];H''-^{n)), it follows that rt G L^{0,T; H'^-^ {fl)) and 

T 

lkt||^2(o-r-H''-i(a)) - (, +'^\\^tt\\L^o.T-M''-^{n))\\r\\LHa,T;H>'in)) 

<P(||r(0)||fe,||r,(0)|U_i) + <5 sup \\r{t)\\l + CT sup (\\r{t)\\l \\ruml-i) ■ 
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Thus, with L'^ -in-time control, we see that the odd time derivatives of rj verify the estimate 



3 

ll^i^°^llL2(0,T;H3 =-»(a)) 
a=0 



<Mo{S)+S sup E{t) + CTP{ sup E{t)) . 
te[o,T] te[o,T] 

See the interpolation inequality {9.21^ below for further details. 

9.2. Assumptions on a priori bounds on [0,Tk]. For the remainder of this section, we 
assume that we have solutions (77^) G on a time interval [0,Tk], and that for all such 
solutions, the time T„ > is taken sufficiently small so that for t E [0,T„] and ^ e M'^, 

det5(rKi))-i/^<2det.gM-i/^-2, \\J-^AIAI-Sl6lh^^n)<h J 
We further assume that our solutions satisfy the bounds 
lhWlll^3.5(a) <2|e|i.5 + l, 

\\d?vmls-.f2^n) ^ 2||9>(0)||^3-./.(a) + 1 for a = 0, 1, 6 , 
IIPo3f??(i)||^,.5_„/.(o) <2||po9?r,(0)||^,.,_/,(f,) + l for a = 0,1,..., 7, 

||V^9,^»+Vt)f^3-.(o) <2||a^(0)|p^3-„(o) + l for a-0,1,2,3. 



(9.4) 



The right-hand sides appearing in the inequalities (|9.4p shall be denoted by a generic 
constant C in the estimates appearing below. In what follows, we will prove that this can 
be achieved in a time interval independent of n. 

We continue to assume that po is smooth coming from our approximation (|7.ip . 

9.3. Curl Estimates. We take T e (0,r«). 
Proposition 9.6. For all t e (0,r), 

X]l!curia^r7(t)llta + E llPo9'-'curl92'^(i)||2 +^ T \\^p,cml, d'-'d?vis)\\lds 

a=0 1=0 1=0 •'^ 

<Mo + CTP{ sup E{t)) . (9.5) 

t6[0,T] 

Proof. Using the definition of the Lagrangian curl operator curl^ given by (|4.ip . we let curl,, 
act on (|7.2a[ ) to obtain the identity 

(curl„ vt)"" = -KekJ^v'^,s [{poJ~^),i K] A™ . (9.6) 

As described above, in the absence of the artificial viscosity term, the right-hand side is 
identically zero; we will have to make additional estimates to control error terms arising 
from K-right-hand side forcing. 
It follows from (HH) that 

5f (curl,, v)'' = ekj^AtJv\s -KekJ^v'',s Aj [{poJ-^),i Al] „„ A™ . 

Defining the fcth-component of the vector field B{A,Dv) by 

B''{A,Dv)^-ek,^A'y,iA'jV\s 
and defining the fc-component of the vector field F by 

F"" = -K£fc,,«^, [{poJ-'),i Al] ,^ AT , 
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we may write 



curl^w(i) = curlMo + / [B{A{t'),Dv{t')) ^ F{t')]dt' . (9.7) 
Jo 

Computing the gradient of this relation yields 

cur\,^Dv{t) ^ D curXuQ - e.j^DA]v\s+ f [DB{A{t'), Dv{t')) + DF{t')]dt' . 

Jo 

Applying the fundamental theorem of calculus once again, shows that 

curl^ D-qit) ^tcmlDuo + e.ji [ [AtjDT]\s -DAjv\s W 

Jo 



t rt' 

[DB{A{t"), Dv{t")) + DF{t")]dt"dt' , 



Jo 



and finally that 



D cnrlriit) ^ tD curl uq - e.ji J Atj {t')dt' Drf 



+ f[At'pT^\,-DA]v\,]dt'+ f f [DB{A{t"),Dv{t")) + DF{t")]dt"dt' . 
Jo Jo Jo 

(9.8) 

Step 1. Estimate for curljy. To obtain an estimate for || curlr7(i)||3, we let act on 
dl^. With dtA"^ = -Afv'',pA''- and DA"^ = -A^Dri^pA^^, we see that the first three terms 

on the right-hand side of (|9.8|) are bounded by Afo + CTP(supjg[Q j-j E{t)), where we remind 

the reader that Mq = P{E{0)) is a polynomial function of E at time < = 0. Since 

DBk{A, Dv) = ~ekjADv\s A^^ + v\, AtDv',p A^^ + v\, Z?(^Mp], 

the highest-order term arising from the action of Z?^ on DB{A, Dv) is written as 

-ekJ^ f f[D^v\sAtv\j,AP + v\,AtD^v\pAP]dt"dt'. 
Jo Jo 

Both summands in the integrand scale like D'^v Dv A A. The precise structure of this sum- 
mand is not very important; rather, the derivative count is the focus. Integrating by parts 
in time, 

/ / D^vDvAAdt"dt' = ^i i D^'q{DvAA)tdt"dt'+i D^rjDvAAdt' 
Jo Jo Jo Jo Jo 

from which it follows that 

pt pt' 

11/ / D'^B{A{t"),Dv{t"))dt"dt'f<CTP{snvE{t)). 
Jo Jo te[o,T] 

We next estimate the term associated to F. Since 

DF'^ - ~K eur. [Dv-,s A] [(po J"'),; K] ,„ A^ + A]D[{p^J-^),i 4] ,„ A^ 

+ v\s[{poJ-'),iA%nD{A]AT)], 

the highest-order term arising from the action of D^ on DF is written as 

^£kl^ f f [DWs A^ [(PO J~'),i A[] AT -f V\s A'jD^ [{poJ-^),l K] ,ra A^W dt' . 

Jo Jo 
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The first summand in the integrand scales hke D'^vD^{poJ^^)AA, and can be estimated 
by integrating by parts in time in a similar way as for the terms associated to D^B(A, Dv). 

The a "priori more problematic term is the second one as it seems to call for five space 
derivatives on 77 that we do not have at our disposal. We first notice that since 

{poJ^^),i K = p,r or] = [Dp]r o 77, 
this integral is under the form 

kI f D^{Dpori)DvAAdt"dt' . 
Jo Jo 

Integrating by parts in time (in the integral from to t'), 

K D'^{Dpor])DvAAdt"dt' = -K / {DvAA)tD^ Dp{ij) dt'" dt" dt' 

Jo Jo Jq Jo Jq 



K [ DvAAD"^ [ Dp{ri)dt" dt' 
Jo Jo 



We now explain why we have control of four space derivatives of the antiderivative (with 
respect to time) of Dp{rj). By definition of the K-problcm (|7.2a[ ). we have 

vt + 2Dpor] + 2K[Dpori]t ^0, (9.9) 

which implies by integrating ()9.9p in time twice that 

2 / /" D'^{Dpoj])dt"dt' + 2k [ D'^iDpoTjjdt' = -D^T](t) +tD^UQ, (9.10) 
Jo Jo Jo 

where we have used the fact that D'^ri{0) = since 77(0) = e. We can now use our Lemma 
3.21 which first yields, independently of k, 

II / f D\Dpor,)dt"dt'\\l<Mo + CEit), 
Jo Jo 

and then by using (I9.10|) . 

t 

KD'^iDp o r])dt'\\l < Mo + CE{t) . (9.11) 
Thanks to (|9.1ip . we get the estimate 

/ D^Fdt"dt'\\l<CTP{ sup E{t)), 
Jo Jo te[o.T] 

and hence 

sup II curl77(<)||^ < Mo + CTP( sup E{t)) . 
te[o,T] te[o,T] 

Step 2. Estimate for curlt^t- From (|9.6p . 

cm-\vt = -e.j, f At]{t')dt' vls+F . (9.12) 

Since 

2Kdt[Dp o 77] + 2Dp orj = —vt , 

by Lemma [221 we see that 

\\[Dpo,^]{t)\\l<Mo + \\vtml, (9.13) 
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from which it immediately foUows that |jF||2 < Mq + CTP(sup(gjQ rp^E{t)). For later use, 
we note from equation (|7.2ap that together with (|9.13p . we have that 

\\Kdt[Dporj]{t)\\l<Mo + \\vtml (9.14) 

Since the highest-order term in D'^B{A,Dv) is D'^v, we see that ||i?(y4, £>v)||2 < Mq + 
CrP(suptg[o,T]^(^)) so that 

II curlut(t)||^ <Mo + CTP{ sup E{t)) . (9.15) 

ie[0,T] 

Step 3. Estimates for curXvttt and curl 9^ v. By time-differentiating (|9.12p . estimating 
in the same way as Step 2, we find that 

II cm\vttt{t)\\l < Mo + CTP{ sup E{t)) , 

te[o,T] 

and 

II cnrld^ v{t)\\l <Mo + CTP{ sup E{t)) . 

te[o,T] 

Step 4. Estimate for pod^ curl 77. To prove this weighted estimate, we write (|9.7p as 

cm\v{t) = ejkzv\s [ At j {f )dt' + curl ua+ [ [B{A, Dv) + F]{t')dt' , 
Jq Jo 

and integrate in time to find that 

curl ry(t) = t curl Mo + / £jktv\sf At]{t")dt" dt' + [ [ B{A,Dv)(t")dt" dt' 
Jq Jq Jq Jq 

" V ' " v ' 

II I2 

+ 11 F{t")dt"dt' . (9.16) 
Jo Jq 

X;, 

It follows that 

pod^ curl77(t) = tpod^ curlug + PqB'^Ii + pQd'^l2 + Pod'^^a ■ (9-17) 
Notice that by definition, 

\\tpod^ cmluoWl < Mo, 

so we must estimate the (fi)-norm of poB^Ii + pod'^l2 + Pod^'^a, ■ We first estimate pod^X2 ■ 
We write 



Po5%(t)=/ / ekjrAt]pod^v\sdt"dt' + / Skj,Pod^Aty,,dt"dt' +n, 
Jq Jq Jq Jo 

ICi IC2 

where TZ denotes remainder terTTis which are lower-order in the derivative count; in particular 
the terms with the highest derivative count in TZ scale like pa^Dv or pd'^rj, and hence 
satisfy the inequality ||7?.(i)||o < Mq + C'TP(supjg[Q E{t)). We focus on the integral /Ci; 
integrating by parts in time, we find that 

JCi{t)=- f f £kJ^^^iA';jPo^\\sdt"dt' + f ek,^At]po^S\sdt' 

Jo Jo Jo 
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and hence 

||/Ci(t)||o<Afo + CTP( sup E{t)). 

telo.T] 

— —ApV^A'j, we see that /C^ 

to yield the inequahty 



Using the identity dtA^ — —ApV^A'^, we see that /C2(i) can be estimated in the same fashion 



\\pod'^l2{t)L < Mo + CTP{ sup E{t)) . (9.18) 

te[o,T] 

Using the same integration- by-parts argument, we have similarly that 

llpoS^WlIn < A^o + CTP( sup E{t)). (9.19) 

te[o,T] 

It thus remains to estimate pod'^I^ in (|9.16l) . Now 

rt ft' 



Po5% = / / pod^F{t")dt"dt' , 
Jo Jo 

the form 

k[ [ poB^Dv D[Dp{T])] A dt"dt' + K [ [ pod'^D[Dp{T])] A Dv dt"dt' +TZ, 
Jo Jo Jo Jo 



10 Jo 

which can be written under the form 



Ti Ta 

where TZ once again denotes a lower-order remainder term which satisfies ||7?.(t)||g < Mq + 

CTP(sup,g[o,j,j^(i)). 
Since 

2KdtD[Dp o 77] + 2D[Dp o 77] = ~Dvt , 
by Lemma [321 we see that independently of k, 

\\D[Dpo7Mt)\\l<Mo + \\vtml<Mo + CE{t) , 

and that 

\\KdtD[Dpo7^]{t)\\l<Mo + CE{t) . 

Thus, the Sobolev embedding theorem shows that 

\\KdtD[DpoTj]{t)\\l^^^-,^^ <Mo + CE{t) . 

Hence, by using the same integration-by-parts in-time argument that we used to estimate 
the integral /Ci above, we immediately obtain the inequality 

\\Timl<A~Io + CTPi sup E{t)). 

te[a,T] 

In order to estimate the integral T2 , we must rely on the structure of the Euler equations 
(|9.9p once again. Integrating in time twice, we see that 

2 / / Pod^D{Dpor]) + 2K [ pod^D{Dp o tj) = -pod^Dr]{t) + tpod'^DuQ . (9.20) 
Jo Jo Jo 

According to Lemma 13.21 independently of k, 



Jo 



Pod^D{Dpori)\\i^^^ .r,L-) < Mo + C\\pod^Dr]{t)U < Mo + CE{t) , 
and then by using (|9.20p . 

\\k j\od^D{DpoT^)\^^^^^^^,^^,^<Mo + CE{t) . (9.21) 
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Returning to the estimate of 72, we integrate- by-parts in time (with respect to the integral 
from to t') to find that 

T2^ f f nf pad'^D[Dp{ri)]{s)ds[A Dv]t{t") dt"dt' 
Jo Jo Jo 

+ f nf pod'^D[Dp{ri)]{t")dt" A Dv{t') dt' . 
Jo Jo 

The inequahty (|9.2ip then shows that 

\\T2ml<Mo + CTPi sup E{t)), 

te[o,T] 

so that 

||po5%(t)||g < Mo + CTP( sup Eit)), 

te[o,T] 

and with (|9TT8l) and (|9?T9| . we see that (|9Tf| shows that 

\\pod'^cvLTlT^{t)\\l < Mo + CT P{ sup E{t)). (9.22) 

te[o,T] 

Step 5. Estimate for pod^ cxxrlvt- From (j9.12p . 

WpoB^ cuTlvtml < \\e■J^Po^' (^j\t^^(t')dt' vU {t)^\\l + \\pod^F{t)\\l 

< CTP{ sup E{t)) + \\pod''F{t)\\l, 
te[o,T] 

and using (19.41) . 

\\pod'F{t)\\l < C\\pod'Dvit)\\l\\nD[Dpirj)]\\l^^^^ + C\\^pod'D[Dpirj)]\\l + CTP{ sup E{t)) . 

te[o,T] 

Since 

\\KD[Dp{r^{t))]\\l^f^,,.^ < C\\KDp{r^{tml < CMo + CT sup E{t) , 

t€[0,T] 

where we have used (|9.14p and the fundamental theorem of calculus. Once again employing 
the fundamental theorem of calculus, 

\\pQd^Dvml<Mo + CT sup E{t), 

te[o,T] 

and hence 

\\pod^Dv{t)\\l\\D[Dp{r^)]\\l^^^-^<Mo + CTP{ sup E{t)) . 

te[o,T] 

On the other hand, since 

2Kpodtd^D[Dp oTj] + 2pod^D[Dp o ry] = -pod^Dvt , 
by Lemma |3.2[ we see that independently of k, 

\\pod^D[Dpon]{t)\\l < Mo + Wpod^'Dvtml < Mo + CE{t) , 

and, in turn, 

\\Kpod^Ddt[Dpor,]{t)\\l < Mo + \\pod^Dvtml < Mo + CE{t). 

By the fundamental theorem of calculus, we thus see that 

\\npod''D[Dp{7])]\\l<Mo + CTPi sup Eit)) 

te[o,T] 
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which shows that 

IIPoS^ curl II 2 < Mo + CTP( sup E{t)) . 

te[o,T] 

Step 6. Estimates for po^^ curl vttt, po5 curl 9^1), and po curl ^Jt;. By time-differentiating 
()9.12p and estimating as in Step 5, we immediately obtain the inequality 

|1po9' CMr\vut{t)\\l + ||po9curiafw(t)||2 + ||po curl9jw(t)||2 <Mo + CTP{ sup E{t)) . 

te[o,T] 

Step 7. Estimate for -^/npo curl^ B'^'v. From (|9.7p 

-\/kpo curl^ = a/kpo^'' curlwo + v^poe-ij^^^ir ^''^j(^) 
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/" y/^pQd^B{A,Dv)dt' + [ y/^pod'^Fdt' +n{t), 
Jo Jo 



where TZ(t) is a lower-order remainder term satisfying \Tl(t)\'^dt < CTP(supjg[o t] ^(0)- 
We see that 

/ WSiWldt' < tMo , 
Jo 

and since \\pQd'^Dr]{t)\\Q is contained in the energy function E{t), 

I \\S2\\ldt' < CTP{ sup E{t)). 
Jo tefo.Tl 



Jensen's inequality shows that 



I WSiWldt' <CTP{ sup E{t)). 
Jo te\o,T] 



te[o,T\ 

The highest-order terms in 54 can be written under the form 

ct ft 



I Kipod^Dv A D[Dp{ri)] Adt' + [ pod"^ D[D p{t])] A Dv Adt' , 
Jo Jo 



Sib 



with all other term being lower-order and easily estimated. By Jensen's inequality and using 



f \\S4a{t')\\ldt' <Ck f t' f \\y^pad'^Dv{t")\\ldt"dt' <CT sup E{t) . 
Jo Jo Jo te[o,T] 

In order to estimate the term S^,, we use the identity 



2K^pod'^D[Dp{rj)]{t)+2y/^ / poB'^ D[D p{7])]dt' = -^pod^Dv{t) 

Jo 

+ y/^poB^Dua + 2k^poB'^D^po , 
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which follows from differentiating the Euler equations. Taking the (f2)-inner-product of 
this equation with poB^ D[D p{ri)]{t) and integrating in time, we deduce that 

f \\Khod^D[Dp{r,)]{t')\\ldt' < Mo+ sup E{t) , 
Jo te[o,T] 

from which it follows, using Jensen's inequality, that 

f WS^bimldt' <CTPi sup E{t)), 

Jo te[o,T] 



and thus 



/ \\y/Kpocmlr,d'^v{t')\\ldt' <CTP{ sup E{t)) 
Jo te[o,T] 



Step 8. Estimates for y^Po curl,, '■df^'v for I = 1,2, 3, 4. Following the identical 

methodology as we used for Step 7, we obtain the desired inequality 



V / \\,/lipocvx\d^-^d'i^v{t')\\ldt' <CT P{ sup E{t)). 
Jo «e[o,T] 



□ 



9.4. K-independent energy estimates for horizontal and time derivatives. We take 

tg (o,r,). 



9.4.1. The d^-prohlem. 

Proposition 9.7. For 5 > Q and letting the constant Mq depend on 1/5, 



sup (\\^od^vml+\\pod^Dn{t)\\l+ f \\^pod^Dv{s)\\lds 
te[o,T] V Jo 



<Mo + d sup E{t)+CTP{ sup E{t)) . (9.23) 
te[o,T] te[o,T] 



FREE-BOUNDARY 3-D COMPRESSIBLE EULER, EQUATIONS IN PHYSICAL VACUUM 59 

Proof. Letting 9^ act on (|7.2ap . and taking the L^(il)-inner product of this with d^v^ yields 
^ Jvl Jn Jn 



Xo Xi X2 

+ K [ d^dta1{plJ-^),kd^v'dx + K [ a>y{pld%J-^),kd\'dx 

Jn Jn 
^ ^ ^ ^ ^ ^ 

X3 X4 

n 



+ K / d^a'lipldtJ-^),kd^v'dx + K / dta'lipld^J~^),kd\'dx 



I5 Is 

3 



/ a^-'af d\'dx + Kj2ci f d^-'dta1{pldKj-^),k d\'dx 

7^l 7^2 

3 

+ acVq / a^-'af (pga'9tJ-2),fe aVd:E+ / afid^pl J-^),kd*v'dx+n5. (9.24) 



The integrals Xa, a — 1,...,6 denote the highest-order terms, while the integrals TZa, 
a = 1, 5 denote lower-order remainder terms, which throughout the paper will consist of 
integrals which can be shown, via elementary inequalities together with our basic assump- 
tions (|9.4I) . to satisfy the following estimate: 

/ na{t)dt<Ma + 5 sup E{t) + CT P{ SM^ E{t)) . (9.25) 
Jo te[o,T] te[o,T] 

The remainder integral TZ^ is comprised of the lower-order terms that are obtained when at 
most three horizontal derivatives are distributed onto /Oq, and although we do not explicitly 
write this term, we will explain its bound directly after the analysis of the remainder term 
7?.4 below. 

We proceed to systematically estimate each of these integrals, and we begin with the 
lower-order remainder terms. 

Analysis of TZi(t)dt. We integrate by parts with respect to Xk and then with respect 
to the time derivative dt, and use (|5.5p to obtain that 

3 x 

7^l = - Vq / / B^^^a'l pIBKJ^^ d'^v\k dxdt 

1 = 1 Jo "'f2 

3 ^ 3 

= ^ Q / [ po {d^-'atdKr^)^p^dS\k dxdt - ^ Q / pod^~'atd'j'^podS\k dx 
1=1 Jo 1=1 

Notice that when ^ = 3, the integrand in the spacetime integral on the right-hand side 
scales like I [dDrj pod^dtJ"^ + dDv p^o^ J~^\ pod*Dr] where £ denotes an L°°{U) function. 
Since \\podfJ~^{t)\\l is contained in the energy function E{t) and since dDr]{t) € L°°{n), 
the first summand is estimated using an L°°-L^-L^ Holder's inequality, while for the second 
summand, we use that ||po>/~^(i)|l4 is contained in E{t) together with an L'^-L'^-L^ Holder's 
inequality. 
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When ^ = 1, the integrand in the spacetime integral on the right-hand side scales like 
£ [dDr] poB^ati + dDv poB^a'^] p^d'^if^k- Since Wp^d'^ Dvt{t)\\f^ is contained in the energy 
function E{t) and since dD-q e L°°(il), the first summand is estimated using an L°°-LF'-L^ 
Holder's inequality. We write the second summand as 

dDvpoB^a^ Pod^ri\p+dDvpod^a\ po5V,3- 

We estimate 



( [ dDvpoB^a^ PodWfi dxdt 
Jo Jn 

= - [ [ [BDvpoB''a^,p poBW + BDv,p poS'af poBW]dxdt 
Jo Jfl 

<C r {\\BDv{t)\\L.^n)\\poB*a{mo WpoBSmL^n) 
Jo 

+ \\B^Dv{t)\\L.^n)\\poB^VmLHn)\\B'a\\o)dt 

<C f {\\BDv{t)\\HO.^(n)\\poBMmo WpoB^vmi 
Jo 

+ \\B^Dv{t)\\Ho.^n}\\poBSmi \\B^a\\o)dt 

<C [ {\\vit)\\H^..^n)\\poB'Driit)\\l + \\vit)\\H2..^n)\\poB'D7j{t)\\oMt)U 
Jo 

+ \\v{t)\\H^.^nMt)\\l)dt, (9.26) 
where we have used Holder's inequality, followed by the Sobolev embeddings 

H°-^{n) ^ L^{n) and H\n) ^ L^{n) . 
We also rely on the interpolation estimate 



\Mh{o.T;H'>-^n)) < C'dl^^WllshlU) ^ + C\\vt\\LHo,T:,mifi))\\v\\L^o.T-HHn)) 

<Mo + d sup Ut)\\l + CT sup (Mt)\\l + \\vtml). (9.27) 

te[0,T] tG[0,T]^ ^ 

where the last inequality follows from Young's and Jensen's inequalities. Using this together 
with the Cauchy-Schwarz inequality, (|9.26p is bounded by CrP(sup(g[Q i?(i)). Next, 
since (|5.6p shows that each component of af is quadratic in Brj, we see that the same 
analysis shows the spacetime integral of BDv poB^a^ poB'^rf' has the same bound, and so 
we have estimated the case ^ = 1. 

For the case that 1 = 2, the integrand in the spacetime integral on the right-hand side of 
the expression for TZi scales like i B^DrjB^Dv p^B'^Drj, so that an — — Holder's 
inequality, followed by the same analysis as for the case I = 1 provides the same bound as 
for the case ^ = 1. 

To deal with the space integral on the right-hand side of the expression for TZi, the 
integral at time t = is equal to zero since r]{x, 0) = x, whereas the integral evaluated at 
t = T is written, using the fundamental theorem of calculus, as 

-Vq / pQB^-'a^^B',r^pQBS\kdx =-Vq /po / {B^-'a1B'j-%poBS\k{T)dx 
Jn *='^ Jn Jq 

which can be estimated in the identical fashion as the corresponding spacetime integral. As 
such, we have shown that TZi has the claimed bound (|9.25p . 
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Analysis of 7l2{t)dt. Using (|5.5I) . we integrate by parts, to find that 



r n2{t)dt = -Y^ci r [ Kd^-'a'i pIB'j- 

Jo Jo Jfi 



' d v\k dxdt 



= ^Cly/K\\d^ 'a'lpod'j ^|lL2(0,T;L2(n)) ll\/Kpo9'*w\fc |lL2(0,T;L2(f2)) 

1=1 

<CTP{ sup E{t)) , 
te[o,T] 

the last inequality following from the fact that ||-\/kpo5'*£'w|||2(o t l2(o)) contained in 
the energy function and that for I — 1,2,3, d'^~^a'l p^d^ contains at most four space 
derivatives of and is controlled L°°-in-time. 

Analysis of TZ^ {t)dt. This remainder integral is estimated in the same way as 7^2 {t)dt. 
Analysis of 7l4{t)dt. Integration by parts using (|5.5p shows that 

'Ri{t)dt^ [ ( a'lid^pl J-^),kd^v'dxdt 
Jq Jn 

^ f {J-^Jt)tdS\kdxdt~ f B^pI J-^a^dS\kdx 

[ B^pI {J-'^Jt)tBS\kdxdt- [ B'^pI B^dWr]{T)dx 
Jn Jn 

b^pI [ dtiJ-^4)dtBS\kiT)dx, 

In Jo 

so that by the Cauchy-Schwarz inequality and Young's inequality, 

[ n^{t)dt < Mo + S sup E{t) +CTP{ sup E{t)) . 
Jo te[o,T] te[o,T] 

Analysis of TZ^{t)dt. The highest-order term is 

/ af(aVo dJ-^),kB^v'dxdt 
iQ Jn 

which can be estimated directly using the Cauchy-Schwarz inequality to yield 
/ 7^5(^)rf^ < Mo + S sup E{t) +CTP{ sup E{t)) . 

Jo t<£[0,T] t<£[0,T] 

Analysis of the integral Io{t)dt. Integrating lo from to T, we see that 

Mt)dt - i / po\B\{T)fdx - Mo . 
^ Jn 

Analysis of the integral Ii{t)dt. To estimate Ii, we first integrate by parts using 



5.5p . to obtain 

Ii= - / B^4pp-'^B'^v\kdx. 
Jn 
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We then use the formula (15. 3p for horizontally differentiating the cofactor matrix: 
Jn 

where the remainder TZ satisfies (|9.25p . We decompose the highest-order term in Ii as the 
sum of the following two integrals: 

Jn 

Iu = - f (aV,.<)(5V,fcaf)dx. 
Jn 

Since = r^t, Xi^j is an exact derivative modulo an antisymmetric commutation with respect 
to the free indices i and r; namely, 

<9S^. atd^v\k 4 = Kd^v\k a'; + (<9S^s < - dS\s K)dWk 4 . (9.28) 

and 

5V,.<5V,fcaf = i^(aV,.a^feaV,a4) - ^aV,s aV,fc i<4)t, (9.29) 

so the first term on the right-hand side of (|9.28p produces an exact time derivative of a 
positive energy contribution. 

For the second term on the right-hand side of (|9.28p . note the identity 

(9477^, c4 - Od^^k a"; = ^J^e,,kdS\r A] e,mnd\'^,s ■ (9.30) 

We have used the permutation symbol e to encode the anti-symmetry in this relation, and 
the basic fact that the trace of the product of symmetric and antisymmetric matrices is 
equal to zero. 

Recalling our notation [curl,,i^] ' = SijkF'^ ,r Aj for a vector-field F, (I9.30p can be written 

as 

(a\^. at - dS\s K)d^v\k a'; = - J' curl,, BSj ■ curl,, d\ , (9.31) 
which can also be written as an exact derivative in time: 

curl,, BS ■ curl, B\ = ^^1 curl„ - B^Tj^r B*rj\, (^JA^^)t + {A^^ADt . 

(9.32) 

The terms in (|9.29l) and (I9.32p which are not the exact time derivatives are quadratic in 
PoB'^Dt] with coefficients in L°°{[0, T] x 17) and can thus be absorbed into remainder integrals 
TZ satisfying the inequality (|9.25p . Letting 

DriB^rj — B^Drj A (matrix multiplication of B'^Drj with A) , 

we have that 

^la = \j f P0^J-'\Dr,BS\^dx-~ I po" J-^\cnY\BW^dx+n. 

where, once again, the remainder TZ satisfies (|9.25p . 

With the notation div,, F = A-fF^,j, the differentiation formula (|5.ip shows that In, can 
be written as 

2 at Jq 



FREE-BOUNDARY 3-D COMPRESSIBLE EULER, EQUATIONS IN PHYSICAL VACUUM 



63 



It follows that 



Ii = ~ [ Po'(J"'|^#^|'- J"Vurl,9Sp- J"Viv,9\P)dx + 7^ 
z at Jq 

Jn 



2dt 



where we have used the fundamental theorem of calculus for the second equality on the 
term Dj^d^r] as well as the fact that was chosen sufficiently small so that ^ < J{t) < |; 
in particular, we write 

= B^Dt] a = d^Drj Id + B^Dt] [ At{s)ds . 

Jo 

It is thus clear that J^^ p^J^^lDj^B'^ijl^dx differs from p^lDB^ril^dx by TZ. Hence, 
j\i[t)ds >\jjo^ (^l^5\(T)p - J-i| curl, aS(r)|2 - J-i| div, BSm?^ dx 
-Mo+ / n{t)dt. 



Analysis of the integral Jg 22{t)dt. B'^^J^^ = —2J^^B'^J plus lower-order terms, which 
have at most three horizontal derivatives acting on J. For such lower-order terms, we 
integrate by parts with respect to dt, and estimate the resulting integrals in the same 
manner as we estimated the remainder term TZi, and obtain the same bound. 
Thus, 



l2 = 2 / p2j-3^r^4^s^^ a^^d\\kdx + n 

Jn 



d 
di 



n 

' pp-^'alBSf,,. a1dS\kdx- [ pl{J'^ala1)t B^rj%r B^ ,k dx + TZ 
n Jn 



^ I plj-^\dYVr,B^r]\^dx + n 

dt Jq 



T n i-T 

2 T-l| j;„ 54^/r7.M2, 



SO that 

( l2{t)dt^ I pp''^\AiY^d'^ri{T)\^dx- Mq+ I n{t)dt . 
Jq Jn Jo 

Analysis of the integral Tj,{t)dt. This follows closely our analysis of the integral Ii. 
We first integrate by parts, using (|5.5I) . to obtain 

Iz = -n [ B^dta^plJ-^ B\\k dx . 
Jn 

We then use the formula l|5.4l) for horizontally differentiating the cofactor matrix: 

l3 = ii f d\\s Ma'; - a>f ] B\\k dx + n, 

Jn 
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where the remainder TZ satisfies (|9.25p . We decompose the highest-order term in I3 as the 
sum of the following two integrals: 

Ja 

1^1, = — n / po^J^^ I div^ S^wpdx . 
Jn 

Since 

Jn 

= K [ poJ'^d\\sAt.d\\kA';dx-- K [ plJ'^\cm\r,d^v\^dx, 
Ja Jn 

and letting 

Drfd'^v = B'^Dv A (matrix multiplication of B^Dv with A) , 
we thus have that 

/ l3{t)dt = n[ I plJ-^\Dr,B^v\'^ dxdt - K I I plJ-^\cw:\r,B^v\'^dxdt 
Jo Ja Jn Ja Jn 

- 1^ [ [ Pa'^J^^ \div^B'^v\'^dxdt+ [ TZdt . 
Jo Jn Jo 



Analysis of the integral Ii[t)dt. Integrating by parts, and using 

( Ii{t)dt:^2K [ [ plJ-^alB'^v^r a\d^v\kdxdt+ [ TZdt 
Ja Jo Jn Jo 

= 2k I [ po^J^^ \div,jB\fdxdt+ [ TZdt. 
Jo Jn Jo 



Analysis of the integral 2z{t)dt. Integrating by parts, and using the Cauchy-Schwarz 
inequality, we see that 



[ I5{t)dt = -K[ [ pldt.r^B^a^B^v\kdxdt 
Jo Jo Jn 



<Cy/K sup \\dtJ ^||L°=(0)||Po9'*a*^||L2(o,T:L2(n))|lv^Po9''w*,fc ||L2(0.T:L2(f2)) 

te[o,T] 

<CTP{ sup E{t)) , 

te[a,T] 

the last inequality following from the Sobolev embedding theorem and the L°°(0, T) control 
of IIpo^VWIIo/ 

Analysis of the integral 2Q{t)dt. Estimating in the same fashion as for Z5 shows that 



/ I(i{t)dt<CTP{ sup E{t)). 
Jo te[a,T] 
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The sum J2l=o lo ^ait)dt. By considering the sum of all the integrals 2ait)dt for 
a = 0, 6, we obtain the inequality 



sup - 

*e[o,T] ^ 



Po\d'^v{t)\^dx+ / plJ-^\d^DT]{t)\^dx] + K 



"'n 



plJ-^\Dr^d\{t)\^dxdt 



<Mo+S sup E{t) + CTP{ sup E{t)) 
te[o,T] «e[o,T] 

+ sup ( plJ-^\d'^ cm\'n{t)\'^dx + K [ [ plJ-^\cm\,^dMt)\^dxdt. 
te[o,T]Jo Jo Jn 

Using the fundamental theorem of calculus, 



pp-'\D,-,d''v{t)\'dxdt^ K / pp-'d''v\rAld''v\sAldxdt 

Jn 



pid\\k dV,k dxdt + K / p^oiJ^'AlAl - Sl5l]d\\r d\\s dxdt , 

Jn 



Jn 



Jn 



where we are using the Kronecker Delta symbol S^, to denote the components of the identity 
matrix. It thus follows from (19.31) that 



sup - 

te[o,T] ^ 



po\d\{t)fdx+ / plJ-^\d^Dr]{t)\^dx 



Jn 



pl\Dd^v{t)\^dxdt 



<Mo + S sup E{t) + CTP{ sup E(t)) 
te[o,T] te[o,T] 

+ sup / plJ~^\d^cm\r]{t)\'^dx + K 
te[Q,T] Jn 



Jn 



pgj^^lcurl^ d\{t)\^dxdt. 



The curl estimates (|9.5p provide the bound for the last two integrals from which the 
desired result is obtained and the proof of the proposition is completed. □ 

Corollary 9.8 (Estimates for the trace of the tangential components of ri{t)). For a = 1,2, 
and 6 > 0, 

sup \ri°'{t)\l^< Mq + 6 sup E + CTP{ sup E{t)) . 

te[o,T] te[o,T] te[o,T] 

Proof. The weighted embedding estimate p.6p shows that 



d'vml<C / pU\dS\' + \d'Drjf)dx. 



Now 



sup / Pald'^ril^dx = sup 
te[o,T]Jn te[o,T 







2 


1 pI 


f d\dt' 


dx < sup Wy^B'^vWl 


In 


Jq 


te[o,T] 



It follows from Proposition 19 . 71 that 



sup \\d'^f]{t)\\l< Mo + CTPi sup E{t)). 



t6[0,T] 



t<£[0,T] 



According to our curl estimates (|9.5p . sup(g[Q || curl?7||| < Afo + C* ^(supjgjo.T] ^(0): 
from which it follows that 

sup \\d^cm\r]{t)\\l,,^y < Mo + CTP{ sup E{t)) , 



te[o,T] 



te[o,T] 



66 



D. COUTAND AND S. SHKOLLER 



since 9 is a horizontal derivative, and integration by parts with respect to d does not produce 
any boundary contributions. From the tangential trace inequality (|6.2p . we find that 

sup |aVWI-o.5 < + CTP{ sup E{t)) , 
te[o,T] te[o,T] 



from which it follows that 



sup \rfit)\l, <Mo + CTP{ sup Eit)) , 
te[o,T] te[o,T] 



□ 



9.5. The 9f -problem. 

Proposition 9.9. For S > and letting the constant Mq depend on 1/6, 

sup (\\^d!vit)\\l + \\podjDvit)\\l+ f \\V^poDd!v{s)\\lds 
te[o,T] V Jo 

<Mo + 6 sup Eit) + CT P{ snp E{t)) . (9.33) 

tG[0,T] t6[0,T] 

Proof. Letting df act on (|7.2ap . and taking the L^(f7)-inner product of this with 9f yields 



1 d 
2dt 



Po\dfv\'dx+ / dfafipt,J~'),kdydx+ / af{pt,d^J~'),kdydx 
n Jn Jn 



Xq Xi X2 

Jn Jn 
+ K f dfa^ipldtJ~^),kdydx + K f dta'^{pld!j-^),kdfv'dx 



I5 le 

7 



Y,ci j dt'aUpldlJ~%kdfv^dx + KY,ci f dt'dtaUpldlJ-^),kdfv^dx 



1=1 1=1 



Kl 112 



+ « Vq / 9f-'af {pldldtJ~^),k dfv' dx . (9.34) 
1=1 -1^ 



The integrals la, a — 1,...,6 denote the highest-order terms, while the integrals TZa, a = 
1, 2, 3 denote lower-order remainder terms, which we will once again prove satisfy 



[ na{t)dt<Mo + 5 sup E{t)+CTP{ sup E{t)) . (9.35) 
Jo *e[o,T] te[o,T] 



The analysis of the lower-order remainder term TZi{t) differs slightly from the corre- 
sponding remainder term in the 9** energy estimates, so we proceed with the details of this 
analysis. 
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Analysis of 7li{t)dt. Using (|5.5|) . we integrate by parts with respect to Xk and then 
with respect to the time derivative dt to obtain that 

7 T 

ni = -y2ci [ I d^-'a'l pldlJ-^ d^v\k dxdt 

7 T 7 

= E ^' / [ {dVJtdlJ-^)^p^dlv\u dxdt -Y^ci I podt'a'ld'j-^podJv\, dx ^ . 
1=1 ^0 Jn Jn 

Notice that when I = 7, the integrand in the spacetime integral on the right-hand side 
scales like t [Dvt podfDv + Dv podjDv] podjDv where £ denotes an L°°{fl) function. Since 
WpqOJ Dv{t)\\Q is contained in the energy function E{t), Dvt{t) is bounded in L°°(il), and 
since we can write podfDv{t) — pQdfDv{0) + p^dj Dv{t')dt' , the first and second sum- 
mands are both estimated using an L°°-L^-L^ Holder's inequality. 

The case Z = 6 is estimated exactly the same way as the case I — 3 in the proof of 
Proposition 19.71 For the case 1 = 5, the integrand in the spacetime integral scales like 
£[DvuPodf J^^ + DvutPoDvtttt]podJ Dv. Both summands can be estimated using an L^-L^- 
Holder's inequality. The case Z = 4 is treated as the case I = 5. The case Z = 3 is also 
treated in the same way as Z = 5. The case Z = 2 is estimated exactly the same way as the 
case Z = 1 in the proof of Proposition 19. 71 The case Z = 1 is treated in the same way as the 
case Z = 7. 

To deal with the space integral on the right-hand side of the expression for TZi, the 
integral at time t = is bounded by Mq, whereas the integral evaluated at i = T is written, 
using the fundamental theorem of calculus, as 

7 7 

Vq / podt'a1dlJ~^podJv\kdx = Vq / podt'aUO)dlJ~^iO)podJv\k{T)dx 
1=1 •'^ 1=1 •'^ 

7 T 

+ Vq /po / {d^t-'a^d\j-\dt' podlv\k{T)dx. 
j^j^ Jn JQ 

The first integral on the right-hand side is estimated using Young's inequality, and is 
bounded by Mq + Jsup^g^Q E(t), while the second integral can be estimated in the iden- 
tical fashion as the corresponding spacetime integral. As such, we have shown that TZi has 
the claimed bound (|9.35p . 

Analysis of J^^ TZ2{t)dt. Using (|5.5p . we integrate by parts, to find that 

rT 7 



f 'R2{t)dt = -y^ci [ f Kdf-^a'l pldlJ-^dfv\k dxdt 
Jo , . Jo Jn 



1=1 

7 



= ^QV^||9f ^a't pad\J ^||l2(o,T;L2(j2)) \\\/^Podtv\k ||L2(o,T;L2(a)) 
1=1 

<CTP{ sup E{t)), 
te[o,T] 

the last inequality following from the fact that Hv^poC^f ^^^|ji2(o T-L^{n)) IIPo9j£'w||^oo(o T-L^(n)) 
are both contained in the energy function. 

Analysis of TZ3(t)dt. This remainder integral is estimated in the same way as TZ2{t)dt. 
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Analysis of the integral J^^ 2oit)dt. Integrating Xq from to T, we see that 

Mt)dt = \l pa\dtv{T)\^dx - Mo . 

Analysis of the integral Ii(t)dt. To estimate Xi, we again integrate by parts using 
to obtain 



Ti=- I dfa'^plJ-^dfv\kdx. 
Jn 



Using the differentiation identity (I5.4p . the same anti-symmetric commutation that we used 
for the (9^-differentiated problem can be employed once again to yield 

pUdJv^,.,At) {dfv\^A^;) = ll|poD,9jz;Wp - ll|poCurl,9j«(OP 

1 

+ 

and 



Hence, 



and 

i-T 



f I,{t)dt f Po^ {\DdJv{T)\^ - J-i| curl, dJv{T)\^ - J'^l div, dJv(T)\^) dx 
Jo ^ Jn 

-Mo+ I n{t)dt, (9.36) 



where the remainder integral TZ satisfies (|9.35p . 

Analysis of the integral l2it)dt. dfj^^ ~ —2J^^dfj plus lower-order terms, which 
have at most seven time derivatives acting on J. For such lower-order terms, we integrate 
by parts with respect to dt, and estimate the resulting integrals in the same manner as we 
estimated the remainder term TZi, and obtain the same bound. 
Thus, 

l2 = 2 / plJ-\ld^,i^\r a'idfv\kdx + n 
Jn 

^Jtl Pp'^^^t^'^r a'ldjv\kdx- j pl{J~^ala^^)tdlv%rdlv\kdx + n 

[ pp-^\div^dJv\'^dx + TZ 
Jn 



- ^ 
~ dt 

so that 



Mt)dt^ / plJ^^\dW,,dJv{T)\^dx- Mo+ / n{t)dt . 
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Analysis of the integral 23{t)dt. This follows closely our analysis of the integral Ii. 
We first integrate by parts, using (|5.5I) . to obtain 

Iz = -n [ dfa'tplJ-^ dfv\k dx . 
Jn 

We then use the formula (|5.4p for horizontally differentiating the cofactor matrix: 

l3 = n [ Po^J'^ dfv'',s [a-a^ - a>- ] dfv\k dx + U, 
Jn 

where the remainder TZ satisfies (I9.25p . We decompose the highest-order term in as the 
sum of the following two integrals: 

Isf, = — K / po^J^^ \divrj dfv\'^dx . 
Jn 

Since 



and letting 



K I poJ~^d^v\sA'^d^v\kA^dx~ K I pgj-i|curl„at^wpdx, 
Jn Jn 



Djid^v — d^Dv A (matrix multiplication of d^Dv with A) , 
we thus have that 

/ l3{t)dt^Kf I poJ^^\D^dfv\^ dxdt - K ( ( plJ-'^\c\xA,^dfv\^dxdt 
Jo Jo Jn Jo Jn 

-K [ [ po^J~^ \AiY ,^dfv\^dxdt+ ( n{t)dt . 
Jo Jn Jo 



Analysis of the integral Ti{t)dt. Integrating by parts, and using (15. 2p . 

I Ii{t)dt^2K [ [ plJ-^ald^v',r a'^d^v\kdxdt+ [ n{t)dt 
Jo Jo Jn Jo 

= 2k f I po^J"^ \AiY,^d^v\^dxdt+ f n{t)dt . 
Jo Jn Jo 

Analysis of the integral Iz{t)dt. Integrating by parts, and using the Cauchy-Schwarz 
inequality, we see that 

[ l5{t)dt = -K[ [ pldtJ-^d^a^d^v\kdxdt 
Jo Jo Jn 

< C^/k sup \\dtJ^'^\\L-'{n)\\podfai\\L2(o^T:L'^(n))\WKPodtv\k \\L^o,T:L^n)) 

te[o,T] 

<CTP{ sup E{t)) , 
te[o,T] 

the last inequality following from the Sobolev embedding theorem and the L°°(0, T) control 
of||po5M(t)llo. 
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Analysis of the integral 2Q(t)dt. Estimating in the same fashion as for I5 shows that 

Ieit)dt <CTP{ sup E{t)). 

te[o,T] 

The sum X]a=o /o^-^a(i)'^^- By considering the sum of all the integrals 2ait)dt for 
a = 0, 6, we obtain the inequality 



1 

sup 7; 

te[o,T] ^ 



pa\d°v{t)\^dx+ / pf,J-^\dlDv{t)\'dx] +K 



<Mo + S sup E{t) + CTP{ sup E{t)) 
te[o,T] te[o,T] 



Jn 



plJ-^\D,jdfv{t)\^dxdt 



+ sup 

te[o,T] Jn 



pp'^ldj cnT\v{t)\^dx + K / / plJ-^lcnrlr, dfv{t)\^dxdt . 

Jo Jn 



Using the fundamental theorem of calculus, 
plJ-^\D,^d^v{t)\^dxdt = K 

■T f 

pldfv\k dfv\k dxdt + K 



Jn 



Jn 



Jn 



pU-'dfv\rAldy,sAidxdt 



pt[J-'AlAl - SlSl]dfv\r dy,s dxdt , 



Jn 



where we are using the Kronecker Delta symbol S^. to denote the components of the identity 
matrix. It thus follows from (19.31) that 



sup - 

te[o,T] ^ 



pa\dfv{t)\^dx+ / plJ-^\djDv{t)\^dx 



Jn 



pl\Ddfv(t)\^dxdt 



<Mo + S sup E{t) + CTP{ sup E{t)) 
te[o,Tl te[o,T] 

+ sup / plJ~^\dJ cuTlv{t)\'^dx + K 



te[o,T] Jn 



Jn 



pIJ^'^I cur\,^dfv{t)\'^dxdt. 



The curl estimates (j9.5l) provide the bound for the last two integrals from which the 



desired result is obtained and the proof of the proposition is completed. 

Corollary 9.10 (Estimates for djv{t)). 

sup \\djv{t)\\l < Mo + S sup £; + CTP( sup E{t)) 
te[o,T] te[o,T] te[o,T] 

Proof. The weighted embedding estimate p.6p shows that 



\\dJv{t)\\l<C pl{\d'v\^ + \DdJv\')dx. 



Now 



pt,\dlv\'dx<Mo+ / p'o 







2 


[ pI 


J dfvdt' 


dx < Mo 


In 







tG[0,T] 

Thus, Proposition 19.91 shows that 

sup \\djv{t)\\l < Mo + S sup ^ + CTP( sup E(t)) 
te[o,T] *e[o,T] te[o,T] 



mi 



□ 



□ 
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9.6. The 9^9"^, dfd^, and dfd problems. Since we have provided detailed proofs of the 
energy estimates for the two end-point cases of ah space derivatives, the problem, and 
all time derivatives, the df problem, we have covered all of the estimation strategies for all 
possible error terms in the three remaining intermediate problems; meanwhile, the energy 
contributions for the three intermediate are found in the identical fashion as for the and 
df problems. As such we have the additional estimate 

Proposition 9.11. For S > and letting the constant Mq depend on 1/S, for a — 1,2, 

sup Y\\drr{t)\i,-a + iivp^s^-" d^vmi + iipo^^-'^ drDm\i 

+ k/ Wpad^-'' d'^''Dv{t)\\lds] <Mq + 5 sup E{t) + CTP{ sup E{t)) . 
Jo te[o,T] te[o,T] 

9.7. Additional elliptic-type estimates for normal derivatives. Our energy estimates 
provide a priori control of horizontal and time derivatives of rj; it remains to gain a priori 
control of the normal (or vertical) derivatives of rj. This is accomplished via a bootstrapping 
procedure relying on having djv{t) bounded in L^{il). 

Proposition 9.12. For t € [Q,T], dfv{t) e H^{n), podfj-'^it) e H^{n) and 

sup {\\d^v{t)\\l + \\podW^{t)\\l) <Mo + S sup E{t) + CTP{ sup E{t)) . 
te[o,T] tG[0,T] te[o,T] 

Proof. We begin by taking six time-derivatives of (I7.2a[ ) to obtain 

2KdJ[AnPoJ"'),k] + 2df[A^ipoJ-'),k] = -djv\ 

According to Lemma and the bound on ||9jti(i)||Q given by Corollarv l9.10[ 

sup \\df[2A'i:{poJ~^),k]\\l< Mo + S sup E{t)+CTP{ sup E{t)). (9.37) 
te[0,T] i6[0,T] te[o,T] 

For /3 = 1,2, 

2A'^{poJ-^),k = PQa'^J-^,k +2po,k a\j-'^ 

= Po4-r\z +2po,3 a- J"' + Poa^J-^p +2po,/3 af J"' (9.38) 

Letting df act on equation (|9.38|) . we have that 

PoaldU^^3 + 2po,3 aldfj-' = 2d! [A'l{poJ-'),k ] - Po^f (af J^^^ ) - 2po,p 9f (af J"') 



Ji ^2 

5 

{dfaf)[poJ-\3+2po,3J-^]+Y.Cadtafdt''[poJ-^3+2po,3 J-^' 

a=l 



Ji 



Jr, 

Bounds for Ji{t). The inequality (|9.37p establishes the L'^iVl) bound for Ji{t). 
Bounds for J2 [t) . According to Proposition 19. Ill 

sup {\\^d!v{t)\\l + \\podDd!v{t)\\l) <Mo + 5 sup E(t) + CTP{ sup E{t)) , (9.39) 
te[o,T] te[o,T] te[o,T] 

so that with (|9.4I) . we see that the differentiation identity (|5.2p shows that 

sup 11:^2(0110 < A^o + sup E{t) + CTP{ sup E{t)). 
te[o,T] te[o,T] te[o,T] 
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Bounds for Jsit). The estimate for J3{t) follows from the mequality for /3 = 1,2 

<C\\po,0h, 



P0,/3 


< c 




< 


PO 




PO 


2 



where the first inequality follows from the Sobolev embedding theorem, and the second from 
our higher-order Hardy inequality Lemma I5TT] since po,p e H^{fl) n -ffp (i^) for (3 — 1,2. 
Thus, 



\\2po,,dKafj-')ro - \\2podnaU-') 

Po 



^e^^P T-2\ P0,I3 ||2 



< \\2podnaU-')\\l 



Po, 



Po 



L°°(n) 



<Mo + 6 sup E{t) + CTP{ sup E{t)) 
te[o,T] te[o,T] 



(9.40) 



thanks to ()9.39p and (|9.4p . and the fact that \\po\\a is bounded by assumption, from which 
it follows that 

sup \\Jz{t)\\l < Mf) + 5 sup E{t) + CTP{ sup E{t)) . 
te[o,T] te[o,T] te[o,T] 

Bounds for ^^{t). The identity (|5.6p shows that is quadratic in drj, and in particular, only 
depends on horizontal derivatives. From the estimate (|9.39p and the weighted embedding 
p.6p . we may infer that 

q|5„.^j.M|2 



Thus 



sup \\dd^v{t)\\i < Mo + 5 sup E{t) + CTP{ sup . 
te[o,T] te[o,T] te[o,T] 



sup ||j4(i)llo < ^0 + ^ sup E{t) + CTP{ sup ^(i)) . 
te[o,T] te[o,T] t<£[a,T] 



Bounds for J^z{t). Each summand in J^z{t) is a lower-order term, such that the time- 
derivative of each summand is controlled by the energy function E{t); as such, the funda- 
mental theorem of calculus shows that 

sup IIJ5WII0 < ^~^o + ^ sup E{t) + CTP{ sup E{t)). 
t<£[0,T] te[o,T] te[o,T] 

We have therefore shown that for all t E [0,r], 

\\poaldW^3+2po,3aldW%<Mo + S sup E{t)+CTPi sup Eit)) , 



te[o,T] 



te[o,T] 



and our objective is to infer that the i^(r2)-norm of each summand on the right-hand side 
is uniformly bounded on [0,r]. 

To this end, we expand the L^(r2)-norm to obtain the inequality 



\\po\a'\d^J-\,it)U + ma'po,3\d^J-'{t)U+4 / poPo,3 lal^r^.V-^S? J'^g dx 



<Mo + S sup E{t) + CTP{ sup E{t)). 
te[o,T] te[o,T\ 



(9.41) 



For each k > 0, solutions to our degenerate parabolic approximation (j7.2p have sufficient 
regularity to ensure that pQ\^^J^'^\'^,z is integrable. As such, we integrate- by- parts with 
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respect to xs to find that 

4 / poPo,3\a'\'dfj'^dfj-',3 dx 
Jn 

= -2\\\^po,3\dfj-^(t)\\l-2 f po{po,3\a'n3(d!j~Ydx. (9.42) 

Substitution of (lO^ into (igTiT]) yieids 

\\po\a'\d!j-',3 ml + 2\\\a'po,3\d!.r'ml 

<Mo + S sup E{t) + CTP{ sup E{t)) + C [ po\df J-^\^ dx . (9.43) 

Using (|9.3p . we see tfiat \a^\'^ has a strictly positive lower-bound. By the physical vacuum 
condition (jl.Sp . for e > taken sufficiently small, there are constants 0i,62 > such that 
Ipoj3 (2:)! > whenever 1 — e < 0:3 < 1, and poix) > O2 whenever < a: < 1 — e; hence, by 
readjusting the constants on the right-hand side of (|9.43p . we find that 

\\p,dU-\3 {t)\\l + 2\\dU-\t)\\l 

<Mq + S sup E{t) + CTP{ sup E{t)) + C [ poldf J'^]^ dx . (9.44) 

By Proposition Em for ;3 = 1, 2, 

sup ||po9f J"^/3 (Oil <^^o + <5 sup E{t) + CTP{ sup E{t)) , 
te[o,T] te[Q.T] te[o.T] 



and by the fundamental theorem of calculus and Proposition l9.91 

sup \\pod^J-^{t)\\ <Mo + 6 sup E{t) + CTP{ sup E{t)) . 
t6[0,T] te[o,T] te[o,T\ 

These two inequality, combined with (I9.44[) . show that 

\\p,dU'\t)\\i + \\dU-\t)\\i 

<Mo + S sup E{t) + CTP{ sup E{t)) + C [ po\df J'^]^ dx . 

We use Young's inequality and the fundamental theorem of calculus (with respect to t) 
for the last integral to find that for 9 > 



c / podtJ-'dU-'dx<e\\dU-Hmo + ^'^\h^t-^ mo 



n+ 

<e\\d!J-\t)\\l + Ce\\pod!Dv{t)\\l 

<0\\dfj-'^{t)\\l + Mo + 6 sup E{t) + CT P{ sup E{t)) , 

te[o,T] te[o,T] 

where we have used the fact that \\podJ Dv{t)\\Q is contained in the energy function E{t). 
We choose 6-^1 and once again readjust the constants; as a result, we see that on [0, T] 

\\podfj~^it)\\l + \\dfj-^it)\\l <Mo + 5 sup E{t) + CTP{ sup E{t)) . (9.45) 

te[o,T] te[o,T] 

With Jt — ajv^jj, we see that 

4 

ald!v\, = dU- v\, d!c4 - J2 Cad^a^dtV,, (9.46) 
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SO that using (|9.45|) together with the fundamental theorem of calculus for the last two 
terms on the right-hand side of (|9.46|) . we see that 

aid^v\j (t) ^ <Mo + S sup £(1) +CTP{ sup E{t)) , 
" te[o,T] te[o,T] 

from which it follows that 

\\divd^v{t)\\l< Mo + S sup E{t) + CTP{ sup E{t)) . 

te[o,T] te[o,T] 

Proposition 19.61 provides the estimate 

\\cuv\d^v{t)\\l<Mo + 5 sup E{t) + CTP{ sup E{t)) 
te[o,T] te[o,T] 

and Proposition 19. Ill shows that for a = 1, 2, 

\d!v°'{t)\l ^<Ma + 6 sup E{t) +CTP{ sup E{t)) . 

te[o,T] te[o,T] 

We thus conclude from Proposition 16.21 that 



□ 



sup \\dfv{t)\\^<Mo + S sup E{t)+CTP{ sup E{t)) . 
te[o,T] ie[o,T] te[o,T] 

Having a good bound for dfv{t) in H^{Q) we proceed with our bootstrapping. 

Proposition 9.13. For t e [Q,T], vttt{t) e H'^{n), padfj-'^{t) e H'^{V) and 

sup {\\vttt{t)\\l + \\p^dtJ-^{t)\\l) <Mo + d sup E{t)+CTPi sup E{t)) . 
te[o,T] te[o,T] te[o,T] 

Proof. We take four time-derivatives of (|7.2a[ ) to obtain 

Kd![A'l{poJ''),k] +2dt[A'^{poJ-'),k]=-d!v\ 

According to Lemma [3^ and the bound on ||c}fw(t)||^ given by Proposition l9.121 

sup \\df[2A'^{poJ-^),k]\\l< Mo + S sup E{t)+CTP{ sup E{t)). (9.47) 

te[o,T] te[o,T] te[o,T] 

For ^ = 1,2, 

2Af (po J"^),fe = Poa-J~^,k +2po,k ap-'^ 

- Pofl? J"',3 +2po,3 (4J-^ + poa^,J-\p +2po,/3 af J"' (9.48) 
Letting ^^ act on equation (I9.48|) . we have that 
poafa,V-2,3 + 2po,3afa,V-2 = 2dt[A^dP^J~').k] ~ podfia^ J-^p) ^2po,^ df{a^ J'^) 

V ' " V ' " V ' 

Ji J72 JTs 

3 

- {dfaf)[poJ~^3+2po,3J-^]+y2cad^afdt-^[poJ-^,+2po,3J-^] . (9.49) 



In order to estimate dfj^^{t) in H^{Vl)^ we first estimate horizontal derivatives of 
dfj^^{t) in i^(ri). As such, we consider for a = 1,2, 

5 

poa^a,V-2,3a +2po,3 a^dfj-^,^ = ^ J(,„ -(poa- ),« ^,^-2,3 -2(po,3 a- ),« J-^ (9 50) 
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Bounds for J7i,q. The estimate (|9.47p shows that 

||J'i,a(t)|lo<^^o + '5 sup Eit) + CTP{ sup E{t)). 

te[o,T] te[o,T] 

Bounds for J2,a- Proposition 19. 1 ll provides the estimate 

sup (^\d'^vutml + IIPo9'^"m(i)||§) < Mo + (5 sup E{t) + CTP{ sup E{t)) . (9.51) 
te[o,T] te[o,T] te[o,T] 

We write 

J2,a = podf{a'^,a J^'^ ,l3a ) + PO,a (9t'(af J"^,^ ) . 

Using (|9.4p . for a — 1,2, the highest-order term in padf{a^J^'^,i3a ) satisfies the inequahty 

which has the bound (|9.51l) . and the lower-order terms have the same bound using the 
fundamental theorem of calculus; for example 

\\PoJ^'^,l3adfa'^J~''^,l3a\\l < 1 1 P0«/"^ , a/3 || L6(f^) || Sfaf || £3 (q) 

< C\\poJ-^,^fi\\i\\dta'^\\o.5 < Mo , 

where we have used Holder's inequality, the Sobolev embedding theorem, and (|9.4p for the 
final inequality. On the other hand, the term pQ,adf{a^J~'^,p ) is estimated in the same 
manner as (19.401) . which shows that 

\\J2,c.{t)\\l<Ma + 5 sup E{t) + CTP{ sup E{t)) . 

te[Q,T] te[o,T] 

Bounds for Jz,a- Using the fact that ||9tJ73,a||o can be bounded by the energy function, 
the fundamental theorem of calculus shows that 

||J3,a Wllo < ^^o + 'J sup E{t) + CTPi sup Eit)). 

te[o,T] te[o,T] 

Bounds for J^,a ■ Again, using the fact that the vector af only contains horizontal derivatives 
of (|9.4p shows that for a — 1, 2, 

II {dta\ [poJ-\^ +2po,3 J"']),a llo < C\\d\ttt\\l + Mo 

<A~Io + S sup E(t) +CTP{ sup E(t)) , 
te[o,T] te[o,T] 

the last inequality following from (j9.5ip . and thus 

\\Ji,c.ml<Mo + S sup E{t) + CTP{ sup E{t)). 

te[o,T] te[o,T] 

Bounds for J^^,a- These are lower-order terms, estimated with the fundamental theorem of 
calculus and (|9.4p . yielding 

\\J5.<.{t)\\l<Mo + 6 sup E{t) + CTP{ sup E{t)) . 

te[o,T] te[o,T] 

Bounds for — (poof ),q dfj^^,3 — 2(po,3 of ),« dfj'^. The bounds for these terms follows in 
the same fashion as for J^2,a and show that 

\\ipoa^),^dtJ-^,3+2ipo,3a^,),a.dtJ-^\\l<Mo+S sup E{t) + CTP{ sup E{t)) . 

te[o,T] te[o,T] 
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We have hence bounded the L^(ri)-norm of the right-hand side of (|9.50|) by Mo+Ssup^^^^ rp-^E{t)- 
CTP(supjg[o,T]^(^))- Using the same integration-by-parts argument just given above in 
the proof of Proposition [9321 we conclude that for a = 1,2, 

sup {\\dtJ~^,c.ml + \\podtJ'^»ml)<Mo + S sup E{t) + CTP{ sup E{t)) . (9.52) 
te[o,T] te[o,T] te[o,T] 

From the inequahty ()9.52p . we may infer that for a =^ 1,2, 

sup \\d\vvtu.c.{t)\\l<Mo + 5 snp E{t) + CTP{ snp E{t)), (9.53) 
te[o,T] te[o,T] te[o,T] 



and according to Proposition [9]6l for a = 1, 2, 

sup II cm\vut,a {t)\\l < Mo + S sup E{t) + CTP{ sup E{t)) . (9.54) 
te[o,T] te[o,T] te[o,T] 

The boundary regularity of Vttt,on = 1, 2, follows from Proposition 19.11} 

sup \vttt,n{t)\l,5< Mq + 6 sup E{t)+CTP{ sup E{t)). (9.55) 

t€[Q,T] t£[0,T] t£[0,T] 

Thus, the inequalities (|g35)) . and (^35)) together with dm and (|g3^ show that 

sup {\\vttt,aml + \\podtJ-\aml) < Mq + S sup E{t)+CTP{ sup E{t)). (9.56) 
te[o,r] te[o,T] te[o,T] 

In order to estimate \\d^J^'^,3, (i)||o, we next differentiate (19.49^ in the vertical direction 
X3 to obtain 



a^S^ J"2^33 +3po,3 a?a,V-2,3 = ^ Ji,3 -poa?,3 -2(po,3 a- ),3 S^V^^ _ (g 57) 



Following our estimates for the horizontal derivatives, the newly acquired inequality (|9.56p 
together with Propositions 19 . 1 ll and 19. l21 show that the right-hand side of (|9.57p is bounded 

in L^i^l) by ^h + 5sup^^^o.T\E{t)+CTP{sup^^^o.T\E{t))■ 
It follows that for k = 1,2,3, 

||poa?afJ-^fc3+3po,3a■9,V-^fc||2 < a1o + <5 sup E{t)+CTP{ sup E{t)) . 

te[o,T] te[o,T] 

Note that the coefficient in front of po,-i a^dfj~^,k has changed from 2 to 3, but the identical 
integration- by-parts argument that we used in the proof of Proposition 19.121 is once again 
employed and shows that 

\\podfJ-\t)\\l + Wdfj-^ml <Mo + S sup E{t) + CTP{ sup E{t)) . 

te[o,T] te[o,T] 

We can thus infer that 

\\AWvttt{t)\\l<Mo + 5 sup E{t) + CTP{ sup E{t)) . 

te[o,T] te[o,T] 

According to Proposition l9.61 || curlu4ti(f)||^ < Mo + 6supifz[Q^T]E{t) + C T P{supf^^o,T]E{t)) 
and with the bound on vf^^ given by Proposition 19 . 1 ll Proposition |6]2] provides the estimate 

\\vttt{t)\\l<Mo + 5 sup E{t)+CTP{ sup E{t)) . 

te[o,T] te[o,T] 

□ 
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Proposition 9.14. For t e [0,T], vt{t) G H^{n), podfJ-'^{t) e H^{n) and 

sup {\\vtml + \\Pod^J^^ml) <Mo + S sup Eit) + CTP{ sup ^(t)) . 
te[o,T] te[o,T] te[o,T] 

Proof. We take two time-derivatives of (|7.2a[ ) to obtain 

Kdf[A'^ipoJ-'),k] +2df[AUpoJ-'U] = -dfv\ 

According to Lemma and the bound on ||9fw(i)||2 given by ProDOsition l9.13[ 

sup \\d^[2A'l{poJ-^),k]\\l< Mo + 6 snp E{t)+CTP{ sup E{t)). (9.58) 
te[o,T] te[o,T] te[o,T] 

Letting df act on equation (19.48^ . we have that 

- {dfaf)[poJ'\3 +2po,3 + CadtafdlpoJ-^s +2po,3 J-^] ■ 

The bound (I9.58P allows us to proceed by using the same argument that we used in the 
proof of Proposition 19.131 and this leads to the desired inequality. □ 

Proposition 9.15. For t e [0,T], ■q{t) e H^{n), pQj~'^{t) e H'^{n) and 

sup {hml + WpaJ'^ml) <Mo+S sup E{t) + CTP{ sup E{t)) . 
tG[0,T] te[o,T] te[o,T] 

Proof. From (I7.2a[ ). we see that 

Kdt[A'l{poJ-^),u] +2[Af(poJ-'),fc] ^-vl. 

According to Lemma |3^ and the bound on ||wt(i)||3 given by Proposition 19. 141 

sup ||2Af (po J"^),fc \\l<Mn + 5 sup E{t) + CTP{ sup E{t)) . (9.59) 
te[o,T] ie[o,T] te[Q,T] 

Since 

we can use the bound (|9.59l) and proceed by using the same argument that we used in the 
proof of Proposition 19.131 to conclude the proof. □ 

10. Proof of Theorem 11.11 (The Main Result) 

10.1. Time of existence and bounds independent of k and existence of solutions 

to pTg]) . Combining the estimates from Propositions [HI O MM ICTI I^T^ IHT^ [^Til 
19.151 and Corollarv 19 . 101 we obtain the following inequality on (OjTk): 

sup E{t) <Mo + S sup E{t) + CTP{ sup E{t)) . 
tG[0,T] te[o,T] te[o,T] 

By choosing S sufficiently small, we have that 

sup E{t) <Mo + CTP{ sup E{t)) . 
te[o,T] te[o.T] 

Using our continuation argument, presented in Section 9 of this provides us with a 
time of existence Ti independent of k and an estimate on (0,ri) independent of k of the 
type: 

sup E{t) < 2Mq , 
te[o,Ti] 
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as long as the conditions (|9.3p and ()9.4p of Subsection [92] hold. These conditions can now be 
verified by using the fundamental theorem of calculus and further shrinking the time-interval 
if necessary. For example, since 

hWll3.5 < 2||e||3.5 + 2 / \\v{t')\\3,5dt' , 

Jo 

we see that for t taken sufficiently small, |l»7(i)ll3 5 ^ 2|ri| + 1. The other conditions in 
Subsection 19.21 are satisfied with similar arguments. This leads us to a time of existence 
T2 > independent of k for which we have the estimate on (0,T2) 

sup E{t) < 2Mo ■ (10.1) 

te[o,T2] 

In particular, our sequence of solutions {77''}k>o to our approximate K-problem (|7.2p satisfy 
the K-independent bound p0.1|) on the K-independent timc-intcrval (0,r2)- 

10.2. The limit as k — 7> 0. By the ^-independent estimate (jlO.ip . standard compactness 
arguments provide the existence of a strongly convergent subsequences for e > 

vf ~^vt inL2((o,r2);i/2(f])). 
Consider the variational form of (|7.2al) : for ah ip e L^{0,T2; H''^{Q)), 



Po{v'^')W'dx^ I pl[r')-\a^'tv\kdx~K I pldt[{J^')-\a^'t]v\kdx 



dt^O. 





The strong convergence of the sequences {t]'^ , ) show that the limit {ri,vt) satisfies 

Povlip^'dx- / Pq J^^af dxdt^O 
Q '-Jn Jn 

which shows that i] is a solution to ()1.9p on the k- independent time interval (0,T2). A 
standard argument shows that v{0) = uq and 77(0) — e. 

10.3. Uniqueness of solutions to the compressible Euler equations ()1.9p . Suppose 
that (77, v) and {fj, v) are both solutions of (|1.9p with the same initial data that satisfies the 
estimate (|1.13p . 
Let 

6v — V — V , 6ri — T] — f] Sa — a — a , SJ^^ — J^^ — , etc. 
Then Sv satisfies 

poSvl + Sa'liplJ-^),k+d'^ip''SJ-^),k = in (0,T] X r!, 

6v^Qon{t = 0}xn. 

Consider the energy function 

4 3 

a=0 a=0 

+ J2 \\pod!^SJ-\t)\\l_, + \\podjDSv{t)\\l + \\pod!dvml ■ 

a=0 
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Given the transport-type structure for the curl of Sr] and its space and time derivatives, 
together with the assumed smoothness of rj and fj, we can proceed in the same fashion as 
our estimates in Section [HI and using that Sv(0) = 0, we obtain 

sup £{t) <CTP{ sup £{t)) , 

tG[0,T] t6[0,T] 

which shows that 5v{t) = on [0,T]. 

10.4. Additional estimates for curl^ v. We have constructed solutions for which 

4 

remains bounded on (0,T). 

The curl estimates required us to assume that the initial velocity field uq satisfies 

II curlwolla < oo and Upo^* curl uqUo < oo . 

Well-posedness requires that the dynamics maintain this regularity for t G (0,T), and we 
show that the Lagrangian curl does indeed maintain the continuity-in-time. 

Corollary 10.1. 

sup (II curl,,w(t)||^ + ||po3'*curl^u(<)|l^) < Mo + ^ sup E{t)+CTP{ sup E{t)) . 

te[0,T] t<£[0,T] t<£[0,T] 

Proof. Letting act on the identity (|9.7p for curl^ v, we see that the highest-order term 
scales like 

curl Wo + / D'^vDvAAdt'. 
Jo 

We integrate by parts to see that the highest-order contribution to D'^ curl^ v{t) can be 
written as 

curl Mo - / D'^ij[DvAA]tdt' + D'^Tj{t)Dv{t)A{t)A{t), 
Jo 

which, according to Proposition 19. l51 has L^{fl)-norin bounded by 

Mo + 5 sup E{t) + CTP{ sup E{t)) , 

te[o,T] te[o,T] 

after readjusting the constants; thus, the inequality for the iJ^(ri)-norm of curl^u(t) is 
proved 

The same type of analysis works for the weighted estimate. After integration by parts in 
time, the highest-order term in the expression for poB^ cml,jv{t) scales like 

po^" curl uo- [ pod'^Dri [Dv A A]tdt' + pod^D7]{t) Dv{t) A{t) A(t) . 
Jo 

Hence, the inequality (19.231) shows that the weighted estimate holds as well. □ 

By taking 6 > and T > small enough, we see that sup^gjo t] ^{^) < ^'^o, where E{t) 
is defined in ([TI^ . 
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10.5. Optimal regularity for initial data. For the purposes of constructing solutions to 
our degenerate parabolic K-problem ()7.2|) . in Section 17.11 we smoothed our initial data so 
that both our initial velocity field Uq is smooth, and and our initial density function pg is 
smooth, positive in the interior, and vanishing on the boundary T with the physical vacuum 
condition (jl.5|) . 

Our a priori estimates then allows us to pass to the limit lim^_>.o Uq = uq and limfl_j.o Po = 
Pq. By construction, po G H*{fl), satisfies po > in fl, and the physical vacuum condition 
(jl.Sp near the boundary F. Similarly, the initial velocity field need only satisfy E{Q) < oo. 

11. The case of general 7 > 1 
We denote by gq the integer satisfying the inequality 

1< 1 + -ao<2. 
7-1 

The general higher-order energy function is given by 
4 4 

a=0 a=0 

+ E llPoSf J-^WIIt, + II curl,t;(0||^ + \\pod'cm\,vml 
+ f2\\^o'^^-''dJ+^-'-^Dvit)\\l, 

a=Q 

and we set = P{E^{0)). 

Notice the last sum in Ej appears whenever 7 < 2, and the number of time-differentiated 
problems increases as 7 approaches 1. We explain this last summation of norms in E-y with 
a particular example; namely, consider the case that 7 = f ■ Then, po ^ <fi near F, ao = 1, 
and the last summation is written as 

"■0 

EllVd ^■^-'''dl+'^^^-'^Dv{t,-)\\l = \\Sd^,Dv{t)\\l + \\d^dlDv{t)\\l 
which is equivalent to 

/ pl\d!Dv{t)\^dx+ [ \pl\djDv{t)\^dx. (11.1) 
The Euler equations with 7 = f are written as 

po«l + af(p|j-i),fe = 0. (11.2) 

Energy estimates on the ninth time-differentiated problem produce the first integral in 
(jll.ip , while the second integral is obtained using our elliptic- type estimates on the seventh 
time-differentiated version of (|11.2p . (Notice that the value of 7 does not play a role in our 
elliptic-type estimates.) Having control on the two integral in (|ll.ip then shows that we are 
back in the situation for the case that 7 > 2; namely, we see that djv{t) is even better than 
L^(r2), which allows us to proceed as before. In particular, for 7 < 2 the power on po in 
the first integral in (|11.1|) is greater than one, and by weighted embedding estimates, this 
means that the embedding occurs into a less regular Sobolev space; this accounts for the 
need to have more time-differentiated problems when 7 < 2. 
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Using this energy function, the same methodology as we used for the case 7 = 2, shows 
that supjg[o T] ^li^) remains bounded for T > taken sufRciently smaU. 

Theorem 11.1 (Existence and uniqueness for the case 7 > 1). Suppose that po £ H'^{il), 
Po{x) > for X (z fl, and po satisfies il.5\) . Furthermore, suppose that Uq is given such that 
Mq < 00. Then there exists a solution to il.9\) (and hence on [0,T] for T > taken 

sufficiently small, such that 

sup E-f{t) < 2Mo. 
te[o,T] 

Moreover if the initial data satisfies 

a=0 

ao 

+ II curl, viO)\\l + WpoB' curl, «(0)||2 + ^ \\^'+^-''d^+^''-^ Dv{t)\\l < 00 , 

a=0 

then the solution is unique. 
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